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INTRODUCTION 


A™ DECADES ago the most enthusiastic admirer of 
mathematical analysis would not have expected 
that practical engineers engaged in the design of air- 
craft would have so much use for the mathematical 
theory of fluid motion as is the case in modern aeronauti- 
calengineering. The various applications of the theory 
are based almost entirely on the hydrodynamics of in- 
compressible perfect fluids, a discipline that half a cen- 
tury ago was considered as an interesting field of pure 
science having very little to do with the motion which 
actually takes place in a real fluid. Recently, interest 
has been centered on another branch of fluid mechanics, 
namely, on the mechanics of compressible fluids. The 
aeronautical engineer is pounding hard on the closed 
door leading into the field of supersonic motion. He 
realizes that the price that has to be paid for further 
increase of speed becomes higher and higher if he nears 
this frontier. However, he wonders whether the mathe- 
matical theory could not be used as a guide for avoiding 
a premature drop of aerodynamic efficiency. The pres- 
ent paper has the objective of reviewing the present 
status of the theory of compressible fluids from the 
practical standpoint of its usefulness for interpretation 
of experimental research and guidance in design. 


IRROTATIONAL MOTION OF IDEAL FLUIDS 


The problem of irrotational motion of ideal incom- 
pressible fluids is more or less a problem of pure ge- 
ometry. Every velocity field that satisfies the condi- 
tions of no vorticity (zero rotation) and the condition 
of continuity (zero divergence) represents a possible, 
1.e., dynamically correct, flow pattern of an ideal in- 
compressible fluid. For example, in the case of two- 
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dimensional motion, any network of two families of 
orthogonal curves represents a network of streamlines 
and equipotential lines of a dynamically correct flow 
pattern provided it satisfies one condition: that all 
its elements are of square shape. The flow is com- 
pletely determined by the rule that the amount of fluid 
passing between two neighboring streamlines is con- 
stant and the ‘‘potential difference’ between two suc- 
cessive equipotential lines is a constant quantity. If 
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the fluid volume passing through an elementary chan- 
nel of width An (Fig. 1) is denoted by Ay, obviously 
Ay = vAn, and if the potential difference between two 
successive cross-sections at the distance As is denoted 
by Ag, obviously Ag.= vAs. If one chooses for con- 
venience Ay = Ay, it is seen that An = As; 7.e., every 
element An As of the network has a square shape. If 
the velocity increases or decreases along a streamline, 
the elements of the network shrink or expand in such a 
way that both As and Am are inversely proportional to 
the velocity. Hence the elements of the network 
always keep their square shapes. 

The pressure does not enter the problem directly, 
because in the case of incompressibility the pressure is 
not related to the physical state of the fluid. This can 
be illustrated, for example, by the fact that provided 
the fluid has no free surface (7.e., a surface of given con- 
stant pressure) one can superpose or remove a gravity 
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(b} 


¢, ¥—networks for flow around a circular cylinder. 
(a) Incompressible fluid. 
(b) Compressible fluid with WM, = 0.400. 


Fic. 2. 


field on the fluid in motion without changing the flow 
pattern or the velocity distribution. The gravity field 
produces a pressure field which is simply superposed 
upon the pressures that keep the fluid in dynamic 
equilibrium. This is correct not only in the case of 
gravity but also in the case of any field of conservative 
forces. 

In dealing with compressible fluids two new aspects 
have to be taken into account: 

(a) The continuity condition for the channels be- 
tween neighboring streamlines involves the variable 
density; namely, one has to write: 


Ay = puvAn (1) 


where Ay is now the fluid mass passing through the 
channel (Fig. 1). The potential difference is equal to 
Ag = vAs as in the case of incompressible fluid, and 
therefore in a network drawn with a constant ratio 
between Ag and Ay the ratio As/ An is ‘variable and 
proportional to the density p. Consequently the ele- 
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ments of the network are elongated in the flow direc- 
tion if the density increases, and expanded normal to 
the flow direction if the density decreases. This is 
clearly seen in Fig. 2 where the ¢, y networks for the 
flow of an incompressible and a compressible fluid 
around a circular cylinder are plotted. 

(b) The flow problem cannot be solved by means of 
geometrical considerations alone. The density is con- 
nected with the pressure and temperature by the equa- 
tion of state. If heat conduction is neglected and it is 
assumed that all particles have started from the same 
state at rest, the density and pressure are connected by 
the “isentropic relation’”’ 


p = const. p” (2) 
where y is the ratio of the two specific heats (for con- 
stant pressure and constant volume) of the fluid. ‘This 
relation, combined with Bernoulli's equation 

vdv + (1/p)dp = 0 (3) 
furnishes a relation between density and velocity. This 


relation has to be introduced into the geometrical con- 
ditions of irrotationality and continuity in order to 
solve the flow problem. 


THERMODYNAMIC RELATIONS 


If Po, po denote the pressure and density of the fluid at 
rest, Eq. (2) becomes 


bP = (po/po")p™ (4) 
On the other hand, Eq. (3) can be written in the form 
vdv + (1/p)(dp/dp)dp = 0 (5) 


The quantity dp/dp is equal to the square of the velocity 
of sound; 1.e., the velocity of propagation of infinitesi- 
mally small pressure changes. With the notations 
dp/dp = a* and (dp/dp)o = y(po/po) = ao” (the velocity 
of sound at rest) the following relations are easily ob- 
tained: 


[(y — 1)/2]v? (6) 
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These relations are well known and are stated here only 
for the reader’s convenience. From (6) and (10) it 
follows that the so-called critical velocity, 7.e., the value 
of v which is equal to the velocity of sound correspond- 
ing to the same pressure and density, and the pressure 
corresponding to this velocity are given by 


Or = V2/(y +1) a0 Per = 2/(y #1)/" po (11) 


For air y = 1.405; therefore v,, = 0.913 ao and p,, = 
0.527 po. 

For moderate values of v/a the relations (8) can be 
approximated by the formulas 


b = po — (pw/2) [1 — (v°/4ao’)] 
and 
b = po — (pv®/2)[1 + (v*/4a*)] 


The last equation can be used, for example, to esti- 


mate the error made by using for the dynamic pressure 


Po — p the value pv?/2 as in the case of incompressible 
fluids. 


MECHANICAL AND ELECTRICAL ANALOGIES 


(a) Consider the flow of a shallow layer of a heavy, 
ideal and incompressible fluid. The fluid is bounded by 
a horizontal plane bottom surface, by vertical walls, 
and by an upper free surface of constant pressure. The 
flow of water in a broad channel or in a river is a good 
example of such a flow, provided the frictional resist- 
ance is neglected or it is assumed that the bottom is 
slightly inclined and the component of gravity in the 
direction of the flow compensates for the frictional re- 
sistance. (This latter assumption can only be realized 
as far as the average resistance is concerned.) If the 
variable depth of the fluid layer is denoted by / and the 
velocity averaged over a vertical line by v, the continu- 
ity equation for stationary flow can be written in the 
form 


hvAn = Ay = const. (12) 
where An is the distance between neighboring stream- 
lines. On the other hand, Bernoulli's equation requires 
that 

(v?/2¢g) + h = const. 

Therefore 
vdv + gdh = 0 (13) 


Denote the quantity gh?/2 by P; then Eq. (13) can be 
written in the form 


udv + (1/h)dP = vdv + (1/h)(dP/dh)dh = 0 (14) 


It is seen that the two-dimensional flow of a layer of 
a compressible fluid of constant depth can be compared 
with the flow of a layer of variable depth of an incom- 
pressible fluid. The Eqs. (12) and (14) are identical 
with the Eqs. (1) and (3) of a compressible fluid if the 


depth h/ is replaced by the density p and the quantity P 
by the pressure ». Asa matter of fact, P is equal to the 
resultant of the pressure forces acting over a section of 
unit width and depth / for an incompressible fluid of 
unit density. Furthermore, dP/dh = gh is equal to 
the square of the wave velocity of long waves propagat- 
ing in a shallow channel of depthh. This wave velocity 
takes the place of the velocity of sound in the gas. - The 
so-called hydraulic jump corresponds to the shock wave. 

However, the relation between pressure and density 
of air is different from the relation between P and h. 
For isentropic change of state, = const. p’, whereas P 
= const. h®. It is seen that the analogy would be per- 
fect for y = 2; however, for air y = 1.40 and the high- 
est value of y (for a monatomic gas) is 1.66. Neverthe- 
less, the analogy is interesting for qualitative informa- 
tion on possible flow patterns, especially in the case of 
supersonic flow of a gas, which corresponds to the so- 
called supercritical flow in a channel or river (‘‘torrent’’). 
Many interesting applications have been made by sev- 
eral authors (cf. bibliography at the end of this paper). 
It seems that E. Jouguet was the first to point out the 
analogy. 

(b) It is known that two-dimensional flow problems 
for incompressible fluids can be solved by use of a so- 
called “electrolytic bath” consisting of a layer of liquid 
of appropriate electrolytic conductivity. Assume that 
an alternating current passes through the liquid, and 
denote the amplitude of the alternating voltage by V, 
the intensity of the alternating current by /, and the 
conductivity of the liquid by o. Then the relation 
between J and the derivative of the voltage in the direc- 
tion of the current vector is given by 


I = o(OV/ds) (15) 


On the other hand, from the continuity of the electric 
flux it follows that 


Ih An = const. (16) 


where / is the depth of the liquid layer and An is the 
normal distance between neighboring lines drawn or- 
thogonal to the lines of constant voltage. 

It is seen that the Eqs. (15) and (16) are analogous to 
the fundamental equations for irrotational flow of an in- 
compressible fluid if V/o is replaced by the velocity 
potential and J by the velocity v. If h is variable, Eq. 
(16) is identical with the continuity equation for a gas 
whose density at an arbitrary point is proportional to 
the value of h at the same point. Hence the “‘electro- 
lytic bath” can be used for solution of problems of com- 
pressible flow by means of an iteration method. One 
can first solve the problem for / = const. by measuring, 
for example, the V distribution for certain given bound- 
ary conditions; then p can be calculated for the entire 
field as function of v, and / modified in such a way that 
h = const. p. By carrying out the experiment for the 
modified shape of the setup one obtains a corrected dis- 
tribution of V and the procedure can be repeated again. 
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G. I. Taylor and C. F. Sharman have carried out such 
experiments; the bottom for the liauid layer was made 
from an insulating material that could easily be carved 
out to an arbitrary shape. Of course, if the gradient of 
h is large, the flow of electricity is no longer two- 
dimensional, and therefore the approximation furnished 
by the electric experiment is not very good. This oc- 
curs especially if the velocity approaches the velocity 


of sound. 


ONE-DIMENSIONAL GASDYNAMICS 


Certain flow problems of compressible fluids can be 
treated as one-dimensional problems. For example, 
the flow of a gas through a duct with slightly variable 
cross-section area can be computed with good approxi- 
mation by introducing mean values for velocity, pres- 
sure and density, the average values being taken over 
cross-sections of the duct. This method of treatment is 
similar to the method generally used in hydraulics. 
Therefore, the theory presented in this section can be 
called ‘“‘gasdynamics in hydraulic fashion.’’ Problems 
in aerodynamics usually do not permit the use of such 
simplified methods; nevertheless, the one-dimensional 
case is of interest to the aerodynamicist, since it gives 
an opportunity to recognize the essential differences 
between the mechanics of incompressible and compres- 
sible fluids. 

If v denotes the mean velocity and p the mean den- 
sity of the fluid in an arbitrary cross-section of area S, 
the flux of mass passing through this cross-section is 
equal to Spv. The conservation of matter or the con- 
tinuity of flow requires that in the stationary case this 
quantity be constant for all cross-sections and that in 
the non-stationary case 

S(Op/dt) + (0/0x)(Spv) = 0 (17) 
The coordinate x is taken along the axis of the duct. In 
the case of stationary flow the velocity and density are 
connected by the relation (5) (Bernoulli's equation). 
From the continuity equation Spy = const. it follows 
that 


(dS/S) + (dp/p) + (dv/v) = 0 (18) 


Substituting from (5) the value of (dp/p) = —(vdv/a?), 
one obtains 
(dp/p) + (dv/v) = [1 — (v*/a*)] (dv/2) 


The ratio v/a is called Mach’s number MM. 
Eq. (18) becomes 


dS/S = —[(dp/p) + (dv/v)] = —(1 — M*)(dv/v) 


It is seen that if 17 < 1 the cross-section area of the duct 
decreases when the velocity increases as in the case of an 
incompressible fluid; whereas for / > 1 the cross-sec- 
tion area increases with increasing velocity. The con- 
dition M = 1 involvesdS = 0. For example, in a con- 
verging-diverging nozzle the velocity of sound can occur 
only in the throat section. 


(19) 


Then 


(20) 
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Relation between the cross-sectional area A and 
the flow velocity in a nozzle. 
(a) Incompressible fluid. 
(b) Compressible fluid. 
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Fig. 3 shows the distribution of cross-section area as 
function of the velocity v, first for an incompressible 
fluid (line a), then for a compressible fluid (line 5) ac- 
cording to Eq. (20). Assume that a certain nozzle 
(Fig. 4) is designed for a pressure ratio (po/p,), where p, 
is smaller than the critical value p,, given by Eq. (11). 
In this case the local velocity of sound is reached in the 
throat cross-section. Then there are two possibilities 
of flow with isentropic change of state. The gas might 
expand to the pressure ~,, flowing with supersonic ve- 
locity, or it might be compressed isentropically. In the 
second case the flow is subsonic and the pressure reaches 
the value p,’ at the exit. The velocity and pressure at 
the exit are the same as in the cross-section of the same 
area in the converging portion of the nozzle. The mass 
flow of gas through the nozzle in the two cases is the 
same. If the exit pressure is greater than p,’, 7.e., it is 
between /) and p,’, the amount of flow varies with the 
pressure difference between entrance and exit. If the 
pressure at the exit varies between p, and p,’, the ca- 
pacity is constant. 

However, the question arises: What happens with 
the gas if the pressure is varied between these limits? 
Observation shows that in this case a more or less sud- 
den change of pressure and velocity occurs somewhere. 
Assume that this shock occurs in a certain cross-section 
of the nozzle and that the change of state downstream 
of the shock is isentropic again; then the gas must have 
a higher value of the entropy than in front of the shock. 
The possible pressure distribution for constant values 
of the entropy which are consistent with the given cross- 
section areas and the given capacity can easily be 
plotted (Fig. 4). The upper branches of such pressure 
curves correspond to flows with subsonic, the lower 
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Fic. 4. Pressure distribution along a Laval nozzle. 


branches to flows with supersonic velocity. It can be 
shown that in the case of such a shock the conditions for 
the conservation of matter, momentum and energy can 
be satisfied only if the transition takes place between 
supersonic and subsonic flow. Hence, a jump is physi- 
cally possible only between points on the initial curve and 
points on the upper branch of one of the other curves. 
The location of the shock is determined by the three 
conditions mentioned. If the end pressure is slightly 
lower than p,’, the shock occurs close to the throat; if 
the pressure is lowered, the shock moves toward the 
exit. There is a value of the end pressure p,” that pro- 
duces a shock right at the exit. If the exit pressure is 
smaller than p,” but greater than p,, the shock occurs in 
the free jet in the form of an “‘oblique’’ shock wave. 

It must be noted that so far the fluid has been con- 
sidered compressible but non-viscous. Therefore, the 
results certainly will be modified by viscosity and espe- 
cially by boundary layer effects. When the nozzle acts 
as diffusor, 7.e., when the flow occurs in the direction of 
adverse pressure gradient, separation may occur and 
change the character of the flow and the pressure dis- 
tribution. However, even if such effects are neglected 
for the time being, the following conclusions are of in- 
terest in view of the general multi-dimensional flow 
problems discussed in the following sections. In the 
case of an incompressible fluid, if the pressure at the 
entrance and exit of a nozzle is arbitrarily given, there 
always exists a flow with a continuous pressure and 
velocity distribution which satisfies the conditions of 
continuity and dynamic equilibrium, the amount of 
flow being a function of the given pressure difference. 
In the case of a compressible fluid this is only true if the 
pressure difference is under a certain limit. Beyond a 
certain limiting value of the pressure difference the 
velocity in the nozzle exceeds the velocity of sound and 
a continuous solution exists only for one definite value 
of the pressure difference. It appears that in the case 
of flows that are partially subsonic, partially super- 
sonic, the existence of continuous velocity fields is one 
of the questions of basic importance. 








STATIONARY FLOw AROUND A Bopy 


The main problem of primary interest for aeronautical 
applications is that of the stationary flow around a body. 
If the body is approximately of cylindrical shape, like 
an airplane wing, the flow can be considered as two- 
dimensional. The flow around a nacelle, fuselage or 
cowling offers examples of typical three-dimensional 
problems. In some cases the assumption of axial sym- 
metry simplifies the problem. 

A multi-dimensional flow can be considered as con- 
sisting of elementary channels of infinitesimal cross-sec- 
tion area. For each of these channels the continuity 
equation and Bernoulli's equation must hold. How- 
ever, since the channels are not separated by solid walls, 
the conditions of dynamic equilibrium must be satisfied 
in directions normal to the flow as well as along the 
channel. If the flow at infinite distance from the body 
is parallel and uniform and the effect of viscosity can 
be neglected, these equilibrium conditions are identi- 
cally satisfied by the existence of a velocity potential, 
1.€., by the fact that the flow is irrotational in the whole 
field. 

Using Cartesian coordinates x, x2, x3, and denoting 
the velocity components by 1, v, v3, the continuity 
equation which corresponds to Eq. (17) can be written 
in the form 


Op ae) 
ar + Dy ae, = 0 (21) 


For stationary flow Ip/dt = 0 and therefore 


“dy ix Op 
— és a on 2 
bein depen”? 


For irrotational flow it follows from Bernoulli's equation 
that forz = 1, 2 and 3 


- > 9% ~ 


p 
Ox; dp Ox; =a Ox; me 2 Ox; 


where v? = v7 + v2? + v3”. Substituting the value of 
Op/dx; from Eq. (23) in Eq. (22) the following relation 


is obtained: 
- Ou; , 1 - - Ov; 
SD (24) 
ff Ox; 3 9 0* 42 F Xk 


Introducing the velocity potential y, defined by the 
relations v; = 0¢/0x;, one obtains 


3 3 3 
d*y 1 dy Op O%y 
= = — = 92 
= fd Ox;? a? > De Ox; Ox,_ Ox;OX, (25) 


It is seen that in the case of an incompressible fluid 
(a = o) the right side of the equation vanishes and ¢ is 
determined by the equation of Laplace. If a is finite, 
the equation is linear in the second derivatives of the 
potential, but the coefficients are functions of the ve- 
locity components and the velocity of sound, which it- 


(23) 
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self is a function of the magnitude of the resulting ve- 
locity according to the relation (6): 


a? = ap” — [(y — 1)/2]v? 


Eq. (25) applied to the two-dimensional problem, 
with the notations x; = x, %2 = y, V1 = Vz, Ve = Vy, can 
be written in the form 


(1 = ae) Be _ 9 Pe (1 = ou) oe = 0 
a® / Ox? a? Oxdy a? / Oy? 

It appears that—with the exception of some cases 
which have very simple symmetry conditions—it is 
extremely difficult to obtain exact solutions of the com- 
plicated non-linear differential Eqs. (25) or (27). In 
the following sections methods tor approximate solution 
of these differential equations are discussed. 


(26) 


(27) 


EXPANSION IN POWERS OF THE MACH NUMBER 


Assume that the velocity potential is expanded in a 
power series of powers of the square of the Mach num- 
ber M in the form: 


¢g = Ul + OM? + O,M*+....] 


where U is the velocity of the undisturbed flow, M@ = 

U/a, and ®o, ,,.... are functions of x1, x2, 3. 
Substituting (28) in Eq. (25) one obtains for the 

coefficient of the zero approximation the equation 


A®, = 0 


(28) 


(29) 


i.e., the equation for incompressible fluids. To obtain 
the coefficient of the next approximation ®, the right 
side of Eq. (25) has to be expanded in powers of M’, 





taking into account the expression (26) for a*. Then 
one obtains 
3 3 
OP, OD, O° 
Ab, = b he & — — (30) 
Ox; Ox, Ox;OX, 


t=1 k=1 


In a more physical language, after solving the flow 
problem for incompressible fluids, the terms on the 
right side of Eq. (30) are being computed using this 
solution and are being substituted in the equation. 
Therefore Eq. (30) has the form 


A®, = F(x, Xe, Xs) (31) 


of a non-homogeneous linear differential equation. If 
®, is known, an equation of similar form is obtained for 
®,, and so on for the higher coefficients. 

The problem of solving Eq. (31) for the given bound- 
ary conditions at the surface of the body is identical 
with the problem of finding the flow of an imcompres- 
sible fluid around the body produced by a field of 
continuously distributed sources of the intensity 
F(x, X2, X3) outside of the body. The flow produced by 
a source around bodies of certain simple shapes like, for 
example, circular or elliptical cylinders, sphere or ellip- 
soids is known and therefore in such cases the solution 
of Eq. (30) can be obtained by superposition. How- 
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Fic. 5. The velocity correction factors vy/vs for flow 
around a circular cylinder. 
(I) First approximation. 

(II) Second approximation, y = 1.405. 

(III) Hodograph method, y = —1, distortion un- 
corrected. 

(IV) Hodograph method, y = .—1, distortion 
corrected. 


ever, the procedure is very laborious, especially for the 
higher coefficients ®2, ®;, etc., and if 7 > 0.4 the series 
(28) converges very slowly. 

The method of expansion—which can be also con- 
sidered as an iteration method—was used by Janzen 
and Lord Rayleigh for computation of the flow of a 
compressible fluid around a circular cylinder; then by 
I. Imai, and by K. Tamada and Y. Saito. 

For the two-dimensional case Imai and Aihara gave 
an elegant method for finding a particular integral of 
the non-homogeneous equation for the first approxima- 
tion: 


Pu ay ea 
m= (= Ox? ai 


4 2B 0b) 0°, 


~ Ox dy Wdy 
(2) 
Oy/ Oy? 
It is known that if o(x, y) is a solution of Laplace’s 
equation, ®p is the real part of the complex function 
F = ®) + 7% of the complex variable z = x + 1y. 
Introduce the complex conjugate of F; namely, F = 
®, — iW and the conjugate of z; namely, z = x — 1y. 
Then the right side of Eq. (32) can be written in the 
simple form R(d?F/dz*)(dF/dz)* and it can be shown 
that the real part of the complex quantity (dF/dz) 
J (dF/dz)*dz is a particular integral of Eq. (32). To find 





(32) 











the complete solution one has to add an integral of the 
homogeneous equation A®, = 0 which together with 
the particular integral obtained satisfies the boundary 
conditions. 

It is seen that provided a closed-form solution ot the 
flow problem for incompressible fluids is known, as, for 
example, in the case of a circular or elliptical cylinder 
or a Joukowsky airfoil, the first approximation for a 
compressible fluid can be obtained by direct integration 
and solution of Laplace’s equation for given boundary 
values. 

Fig. 5 refers to the flow around a circular cylinder. 
It gives the correction factor to be applied to the ve- 
locity distribution over a section passing through the 
center of the circle normal to the main flow direction; 
y is the distance from the center, ry the radius of the 
cylinder. The correction factor plotted in the diagram 
is vy4/Uo, Where vy refers to M = 0.4 and w to the solu- 
tion for incompressible fluids. The two curves (I) and 
(II) show the first and second approximations obtained 
by Imai and Tamada by using the method of sources 
(developed by L. Poggi). The value y = 1.405 is used 
for the second approximation; the first approximation 
is independent of the value of y. The rest of the curves 
plotted in the figure will be discussed in the section on 
Applications of the Hodograph Method. 

The following table shows the particular cases treated 
by means of the method of expansion up to the present. 





First Second 
Approximation Approximation 

Circular cylinder Lord Rayleigh, O. Jan- I. Imai, K. Tamada 

zen, L. Poggi and Y. Saito 
Elliptic cylinder S. G. Hooker, C. Kap- 

lan, I. Imai and T. 

Aihara 
Joukowsky airfoil L. Poggi, C. Kaplan 
Sphere K. Tamada K. Tamada 





METHOD OF SMALL PERTURBATIONS 


It is known that the so-called thin airfoil theory has 
been applied with success to airfoils of small thickness. 
This theory is based on the assumption that the velocity 
changes caused in a uniform parallel airstream by the 
presence of the airfoil are small in comparison with the 
velocity of the undisturbed flow. 

The same assumption was applied by Glauert, 
Ackeret and Prandtl to the flow of a compressible fluid 
around airfoils or slender bodies. Denote the velocity 
of the undisturbed flow by V and assume that the 
squares and products of the quantities (v, — V), v2 and 
v3 can be neglected. Then Eq. (24) obtains the form 


a? Ox, (33) 


The local velocity of sound is determined by the rela- 
tion of Eq. (6). 
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It is seen that the difference between a? and a,?, 7.e., 
between the squares of the local velocity of sound and 
the velocity of sound in the undisturbed fluid can be 
neglected in this approximation. With the introduc- 
tion of the velocity potential in the form g¢ = Vx + ¢’ 
Eq. (33) furnishes the following linear differential equa- 
tion for the ‘‘perturbation potential’’ ¢’: 


V2\ d20’ d%y’ d2y’ 

(1 Y) Ox” + Ox»? + x32 ’ 
This equation can be reduced to Laplace’s equation 
provided V < a, by introducing as new independent 
variables & = , & = %*/1 — M;? and & = 
x3\/1 — M,2, where M, = V/a; is the Mach number for 
the undisturbed flow. If (V/a,) > 1, 2@.e., the undis- 
turbed flow is supersonic, Eq. (35) has the character of 
the wave equation. 

From Eq. (35) it is seen that the potential function 
for a flow in an incompressible fluid around a surface 
whose equation is given by x2 = f(x:), x; = g(x) can 
be expressed approximately in the form 


= Va + 9'(%1, V1 — M3, xV/1 — M;?) 
(36) 








(35) 





(x1, X2, X3) 


where ¢’ is a solution of the Laplace equation for the 

variables & = x1, & = x2\/1 — Mi, & = xx/1 — M,?. 

This result can be interpreted in the following way. 
Evidently 


o(éi, f, &) = Vii + 9’ (&, &, &) 


represents a flow of an incompressible fluid in the 
£1, £, &; coordinate system. One has then to investigate 
the shape of the body in this system by comparison of 
the streamlines in the x1, x2, x3 and the §£;, £, £; coordinate 
system. For sake of simplicity take a body with axial 
symmetry and compute the inclination 6 of an arbitrary 
streamline in the x; x2 plane. With the approximation 
used in this section tan B = (1/V)(Og’/Ox2) in the 
X1, %2, X3 system and tan 8’ = (1/V)(Oe’/O&) in the 
&, &, & system. It is seen that 


as os 
Ox, V1 — Mi? 


From Eq. (37) it follows that the flow in the &, &, & 
system represents the flow around a body obtained from 
the original body by expanding its lateral dimension in 
the ratio 1/1 — M,?. 

Therefore, if one has to determine the flow of a com- 
pressible fluid around a slender body, the Mach number 
of the undisturbed flow being equal to Mj, this can be 
accomplished approximately by expanding the dimen- 
sions of the body normal to the main flow in the ratio 
1:\/1 — M,? and solving the problem of the flow of an 
incompressible fluid around the body of the expanded 
shape. Especially it can be shown that with the ap- 


tan B’ _ Oy’ 
tan B ‘i Of 








(37) 
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proximation used in this section the pressures acting on 
corresponding points of the two bodies are equal. 
Namely neglecting the squares and products of per- 
turbance velocity components considered small in com- 
parison with V the pressure at an arbitrary point of a 
compressible or incompressible fluid is equal to 


bp — fi = —aVin — V) = —pVr’ 


where v,’ denotes the perturbation velocity in the direc- 
tion of the main flow and /, refers to the undisturbed 
flow. Since (Og’/Ox2) = (Og’/Ok:), z.e., the v,’ compo- 
nents are equal in the two cases, the pressures are also 
approximately equal. 

Compare now the pressure distribution acting in an 
incompressible fluid on two slender bodies, assuming 
that the shape of the second is obtained by lateral ex- 
pansion in the ratio 1:. If the potential function for 
the flow around the first body is 


g= Vi t+ ¢’ (1, fo, §3) 


the potential function for the flow around the second 
body is approximately given by 
g = Vir + ne'(&s, &2, &5) 


This can be shown by consideration of the inclination of 
the streamlines in the manner indicated above. Then 
it is seen from Eq. (38) that the perturbation velocities 
v,’ and the pressure differences p — ); are also in the 
ratio 1:n. 

This result can be applied to an approximate calcu- 
lation of the pressure coefficient in a compressible fluid 
if the value of the pressure coefficient for the flow of an 
incompressible fluid around the same body is known. 
It was shown that the pressure acting in a compressible 
fluid is equal to the pressure acting in an incompressible 
fluid on a body obtained by expansion in the ratio 
1:\/1 — M,2. The incompressible flow around this 
body can now be compared with the incompressible 
flow around the original body. As it was explained 
above, the pressure coefficients vary proportionally with 
the lateral dimensions; hence the pressure coefficients 
at corresponding points are in the ratio 1:\/1 — M,?. 
Therefore, if the pressure coefficient at an arbitrary 
point of a body in an incompressible fluid is equal to 
Cre, = Ap/(p/2)V*, the pressure coefficient in a com- 
pressible fluid is approximately equal to Coy = 
(1/V/1 = M;*)C»,. 

This relation of course cannot be correct for the 
neighborhood of the stagnation point, where the perturb- 
ance-velocity is of the same order as the undisturbed 
velocity V. Nevertheless, if the influence of such re- 
gions is small, the same rule applies to the total forces 
acting on the body. For example, as was shown by 
Glauert, Ackeret and Prandtl, the lift coefficient of an 
airfoil moving with high velocity can be obtained from 
the lift coefficient obtained for small values of M, by 


the formula : 
Cry = (/V1 — My*)Cz, 


(38) 


(39) 


(40) 
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In the case of an airfoil the ‘‘expansion’’ normal to the 
main flow is essentially equivalent to a change of angle 
of attack, camber and thickness. 

According to Eq. (40) the slope of the lift curve 
dC,/da should also vary in proportion to 1/1/1 — M,?. 
However, the comparison with experimental results does 
not permit one to formulate convincing final conclu- 
sions. In some cases good agreement is found; in 
other cases the rule is at variance with the experimental 
evidence. It is the author’s impression that in most 
cases the Glauert-Prandtl rule underestimates the com- 
pressibility effect on the slope of the lift curve. For 
example, in Fig. 6 the ratio (dC,/da)y/(dCz/da)o is 
plotted as a function of M,. The experimental points 
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Fic. 6. Relation between lift-curve slope dC,/da and 
Mach number M, for N.A.C.A. 4412 (J. Stack, W. F. 
Lindsey, R. E. Littell). 
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are computed from Fig. 9 of N.A.C.A. Technical Re- 
port No. 646 and refer to the N.A.C.A. 4412 section. 
The value for (dC,/da)o is taken from the pressure 
distribution tests at alow Mach number. Fig. 7 shows 
the values of the same ratio for a propeller section of 
6.32 per cent thickness investigated in the high-speed 
tunnel at Guidonia. 

The slope dC,/da reaches a peak value at a Mach 
number between about 0.6 and 0.8; beyond this peak a 
more or less rapid drop takes place. This drop appar- 
ently coincides with the first occurrence of a shock wave. 
This effect will be discussed in a later section. 


THE HopOGRAPH METHOD 


The analytical method presented in this section is re- 
stricted to two-dimensional problems. In the section 
on One-Dimensional Gasdynamics the fundamental dif- 
ferential equation for the two-dimensional irrotational 
motion of a compressible fluid in Cartesian coordinates 
is given. In many problems it is convenient to use the 
are length along streamlines, and equipotential lines as 
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Fic. 7. Relation between lift-curve slope dC,/da and 
Mach number M, for a 6 per cent thick Italian propeller 
section (A. Ferri). 


independent variables; the magnitude of the velocity v 
and the angle of inclination between the velocity vector 
and a fixed direction (e.g., the direction of the undis- 
turbed flow) as unknown functions. Then the velocity 
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distribution is determined by the conditions of irrota- 
tionality and continuity. The rotation of the fluid in 
the case of curved flow is given by (Fig. 8) 


w = (Ov/dn) + (v/R) 


where v is the magnitude of the velocity, is the arc 
length taken along an equipotential line, and R is the 
radius of curvature of the streamline. Evidently 1/R 
= 00/0s, where s is the arc length taken along a stream- 
line. Hence, one obtains as the condition of irrota- 
tionality : 


(41) 


-——— =0 (42) 


The continuity condition requires that for an ele- 
mentary channel between two neighboring streamlines 
pvuAn = const., 2.e., 


1 Op 1 Ov 1 OAn 

oe eape — —— = 43 

p Os v Os An Os : ad 
According to Eq. (19) 

10 . lo », 1 Ov 

a oe ae OR pi 44 

a ve ql OS ts (44) 


345 
On the other hand, it is seen from Fig. 8 that 
1 dan _ 20 ie 
An Os On 
Hence one obtains the equation 
1 Ov rele) 
(1— M*?)-—+—=0 (46) 
vos On 


The two Eqs. (42) and (46) determine the velocity 
distribution in the following sense. Assume one draws 
arbitrarily an orthogonal system of curves g = const., 
¥ = const. Then by integrating Eqs. (42) and (46), 
for example, by a step by step method, one obtains cer- 
tain values for v and 6. Obviously the values of the 
angle © obtained in this way should be consistent with 
the arbitrarily assumed streamlines, 0 indicating the 
direction of the velocity vector and therefore the stream- 
line at an arbitrary point. If this is not the case, one 
has to correct the network and continue the procedure 
which is essentially a trial and error method. The Rus- 
sian mathematician, Chaplygin, noticed that the prob- 
lem can be greatly simplified if a change of variables is 
performed by using v and 0 as independent variables 
and g and y as unknown functions. His method 
amounts essentially to plotting a network ¢ = const., 
y = const. in a plane in which the velocity components 
v, and v, are the Cartesian coordinates or v and 0 the 
polar coordinates. In other words, as one proceeds 
along a streamline or an equipotential line in the “‘physi- 
cal plane’”’ (the plane in which the flow actually takes 
place), one plots (Fig. 9) the velocity vector as to magni- 
tude and direction starting from the origin in a second 
% y 

Y-const 
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HODOGRAPH PLANE 


plane, the ‘‘hodograph plane.’’ The curves described 
by the end points of the velocity vector are considered 
as the curves y = const. and g = const., respectively. 

It is important to obtain the equations for ¢ and y 
as functions of the new independent variables v and 0. 
The first step consists of introducing in Eqs. (42) and 
(46) y and y as independent avariables instead of s and 
n by means of the relations dp = vds and dy = pvdn. 
Then Eqs. (42) and (46) obtain the form 


pov 00 _ 
v Ow de 
1 Ov 00 
~ M5 -— $o— = 47 
(1 sae te 0 (47) 
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Now the independent and dependent variables have to 
be interchanged. This is a well-known operation in 
mathematical analysis. If y is a function of the inde- 
pendent variable x and one chooses to consider x as a 
function of y, obviously from dx/dx = (dx/dy)(dy/dx) 
= 1 follows dx/dy = 1/(dy/dx). Using the obvious 
identities: O¢/O0¢ = Oy/dy = 1, 0¢/dy = OY/de = 0, 
the following analogous relations between the partial 
differential coefficients are obtained: 


ee ee 
oy Ade’ dy = Ado 
ee ee. ee 
%e jjjj.Adve o% Ad (48) 
where A is the so-called ‘‘Jacobian’’; defined by 
_ 262%) d¢dy ro 
ov 00 00 Ov 


The expression A is denoted in general by O(y, y)/- 
O(v, 8) and has to some extent properties analogous to a 
differential quotient. For example, the product of dif- 
ferentials dgdy is equal to (0(y, Y)/O(v, 8))dvdO, just as 
the differential dy is equal to (dy/dx)dx. 

Introducing the expression (48) in Eq. (47) one ob- 
tains: 


poe OW _ 0 
v 00 Ov 
1° . o 


i. es Se. a. 
a oa. oe , (50) 


(51) 





Eliminating ¢ the calculation results in obtaining the 
following linear partial differential equation for y as a 
function of v and 0: 
—_ 2 2 
1 wee. (22) — 
p” oe? p Ov 

It is evident that since v and 0 are the independent 
variables, and p and M are functions of v, this equation 
is a linear differential equation with variable coefficients, 
but the variable coefficients are functions of the inde- 
pendent variables only whereas in the case of Eq. (27) 
the coefficients are functions of the derivative of the 
unknown variable. Therefore, the solution of the Eq. 
(27) can be obtained only by a laborious iteration proc- 
ess, whereas the methods of analysis can be applied 
to Eq. (51) for obtaining families of exact solutions. 
This does not mean that one is able to solve exactly the 
flow problem for given boundary conditions in the physi- 
cal plane, but one obtains exact solutions in the hodo- 
graph plane, which can be transferred to the physical 
plane and furnish certain exact flow patterns of a com- 
pressible fluid. 

Chaplygin and his followers noticed that Eq. (51) is 
reduced to Laplace’s equation if it is assumed that 
(1 — M*)/p? is a constant quantity. Using the general 
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thermodynamic relations, Eqs. (6) and (7), this expres- 
sion is obtained in the form 


1— M 1 . 
— -30-4 “) x 











p* po” 2 ao’ 
~ 1 y?\Gt+n/a- 
ee i (52) 
2 Ao? 
It is seen that 1 — M?/p? = const. if y = —1. For 


this value of y the relation between density and velocity 
is given by Eq. (7) 

p = poll + (v?/ao?)] ~ 
and according to Eq. (6) the square of the velocity of 
sound is equal to 


(53) 


a? = ao*[1 + (v/a*)] = ao? (07/p?) (54) 
Now using the definition of a? 
a® = dp/dp = ao"(po°/p*) (55) 


one obtains by integration the pressure-density relation 


p = const. — (dop9?/p) (56) 


This relation, of course, is not satisfied by any real gas, 
since for all real gases 1 < y < 1.66, however, the con- 
stant of integration in Eq. (53) can be determined in 
such a way that in the p, 1/p diagram (Fig. 10) the 
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Fic. 10. Approximations to true isentropic pressure- 
volume curve. 


straight line defined by Eq. (53) is tangent to the actual 
isentropic curve of the real gas. Chaplygin determined 
the constant so that the straight line approximates the 
isentropic curve near the point corresponding to the 
state of the gasinrest. In this way he and his followers 
solved certain problems by integration of Eq. (51) and 
obtained good approximations if the maximum velocity 
did not exceed about 0.4 times the velocity of sound. 

In order to apply the hodograph method to larger 
velocities one can determine the constant of integration 





COMPRESSIBILITY 


in Eq. (55) so that the isentropic curve of the gas is 
approximated near the point ,, 1/p; which corresponds 
to the undisturbed flow with velocity V. This idea was 
carried out at the suggestion of the present author by 
H. S. Tsien, and will be discussed in the following sec- 
tion. 


APPLICATIONS OF THE HODOGRAPH METHOD 


If the quantity (1 — M*)/p? is assumed constant, 
after division of Eq. (51) by (1 — M?)/p?, this constant 
factor can be included in the bracket and one obtains 


2 oe 
mo + (yrtaaa) #o° 


Assume that w, 6 denote the magnitude and the angle 
of inclination of the velocity vector of an incompressible 
fluid. Then for an arbitrary flow pattern of this in- 
compressible fluid the stream function satisfies the equa- 
tion 


(57) 


0*y/00* + [w(0/ow)]*y = 0 (58) 


Consequently, if ¥(w, @) is a known flow pattern of an 
incompressible fluid and w and v are connected by the 
relation 


V1 — M%(dv/v) = (dw/w) (59) 


the function ¥(w(v), 6) determines the flow of a com- 
pressible fluid characterized by the pressure-density 
ratio (56). The factor 7/1 — M? can be expressed by 
the local velocity v, the velocity V and the Mach num- 
ber M4, of the undisturbed flow. Taking into account 
the relation 


a2 — vy? = a,? — J? 


where d; is the velocity of sound in the undisturbed flow, 
one obtains 
v? vy? 


M? = — = —— = 


— M,*(v/V)? 
ee a°- V+ 


1 — Mi? + M,(0/V)? 





and 





V1 — M2? = V1 — M?/V1 — M2 + M°(0/V)? 


Hence, the relation between the velocities of the in- 
compressible and compressible fluids is given by 





dw/w = V1 — MY/V1 — M2 + M2(0/V)%d0/2) 
(60) 
By integrating and putting w (l’) = W one obtains 
* 4 l+Vi- Me 
Wo VVi— Mi + M20/V + V1 — Me 
or 
v/V = (w/W)(1 — X)/1 — A(w/W)? (61) 
where 
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= M2/(1 + V1 — M,?)? 
The pressure coefficient corresponding to the Mach 
number M, = V/a, is given by 








_ Ap — pi) _ 
Ong = AVE 
, , 
C (62) 
is Mi? Cy, 


Vi — M2+ 





1+ V1 — M; 2 


provided C, is the value of the pressure coefficient for 
M,= 0. Itis seen that Eq. (62) differs from the Glau- 
ert-Prandtl formula by the second term in the denomi- 
nator which itself contains C, again. Hence, if C, is 
negative, the magnitude of the compressibility effect 
predicted by Eq. (62) is greater than that predicted by 
the Glauert-Prandtl formula. The formula (62) can 
be applied to calculation of the influence of compressi- 
bility on the pressure distribution of airfoils, if the pres- 
sure distribution at low Mach numbers is known. In 
this application a slight distortion of the shape of the air- 
foil is disregarded; namely, if ¥(w, 0) is the stream func- 
tion for the flow of an incompressible fluid around a 
given section, so that, for example, the streamline y = 
0 corresponds to the boundary of the section, then the 
streamline ¥y = 0 of the corresponding compressible 
flow in the physical plane will not coincide exactly with 
the boundary of the original section. It appears—at 
least in an example investigated by Tsien and discussed 
later in this section—that the distortion of the section 
is relatively small (with the exception of such thick sec- 
tions as a circle) and the error made in the pressure coef- 
ficient tends to correct the error made by use of the 
approximate pressure-density relation. 
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Fic. 11. Pressure distributions for N.A.C.A. 4412 at 
° 


a= —2°, 


In Fig. 11 the measured pressure distributions of the 
N.A.C.A. 4412 airfoil are plotted for M, = 0.100 and 
M, = 0.606 according to N.A.C.A. Technical Report 
The dotted line is calculated for M, = 0.606 
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by Eq. (62) using for C, the values measured at M, = 
0.100 and a = —2 degrees. 

Figs. 12 (a, b, c) show the variation of pressure coeffi- 
cients at three points of the same airfoil section; the 
measured values are compared with the values computed 
by Eq. (62) and by the Glauert-Prandtl rule. 

Strictly speaking, the theory developed in this section 
deals with a specific kind of compressible fluid for which 
7 = —1. However, one can look at the method also 
from a different point of view. For air, the factor 
(1—M?)/p? in Eq. (51) is variable, but a fair approxima- 
tion might be expected if this variable coefficient is re- 
placed by its mean value. This is probably correct 
unless the local Mach number comes close to unity. 
Replacing (1 — M?*)/p? by its value corresponding to 
the undisturbed fluid, 7.e., by (1 — M,*)/p:”, one obtains 
the following equation for y: 


oy ( P1 2)’ 
= yee = 0 63 
00? p V1 — M,2o . oe 
Then Eq. (59) is replaced by 
d d 
V1 —- Mee = — (64) 


Eq. (64) leads to the following approximate relation 
between v and w: 
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(w/W)* ad (v/ Vg Sr /ate/9 —1) (65) 
where 
y-1 —i/(y¥—-1) a 
b= (1+ : Ms) /V1 — M? 


The pressure coefficients calculated from (65) are shown 
in Fig. 13; they are only slightly different from those 
obtained from Eq. (62) plotted in the same figure by 
dotted lines. For larger values of M, Eq. (65) appar- 
ently gives better agreement with the experiments. 

The physical meaning of the Eqs. (63) and (57) can 
be described in the following way. The stream func- 
tion y[w(v), 8] which is a solution of Eq. (63) furnishes 
a family of approximate streamlines for a compressible 
fluid satisfying the exact relations between pressure, 
density and velocity; whereas the stream function satis- 
fying Eq. (57) furnishes exact streamlines of a compres- 
sible fluid whose physical behavior is only approxi- 
mately correct. The second method has the advantage 
of simplicity in the mathematical expressions. 

In order to have an estimate of the error made by use 
of the approximate pressure-density relation, a calcula- 
tion was made in an extreme case; namely, in the case of 
a circular cylinder. Fig. 5, which was partly discussed 
before in the section on Expansion in Powers on the 
Mach Number, shows the velocity correction factors 
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Fic. 12 (c). Increase of the pressure coefficient C, with 
Mach number M, for N.A.C.A. 4412. 
At 2.5 per cent chord, lower surface, a = —2°. 


vy/v for M, = 0.4 calculated under two different as- 
sumptions. Curve III is obtained without correction 
of the distortion of the original section, Curve IV by 
solving the flow pattern for an originally elliptical cylin- 
der which by the distortion becomes almost exactly 
circular (within '/50 of a per cent of the radius). In both 
cases the approximate pressure-density relation (56) is 
employed, which makes an exact solution of the flow 
equations possible. These curves can be compared 
with the second approximation obtained for y = 1.405 
(Curve II). This curve probably lies somewhat below 
the curve which would represent the exact solution for 
y = 1.405. It is seen that the error made by the dis- 
tortion lies in the right direction to balance the error 
made by using y = —1 instead of the correct value of y. 


DRAG AND TRANSITION THROUGH THE VELOCITY OF 
SOUND 


The experimental evidence shows that the compressi- 
bility effect on the drag of airfoil sections can be attrib- 
uted to two different physical processes. At low val- 
ues of the Mach number the drag of most airfoil sections 
consists essentially of skin friction of partly laminar 
and to a greater part turbulent character. The com- 
pressibility effect on the skin friction coefficient of a 
boundary layer with a fixed pressure gradient is very 
small in the range of the Mach numbers considered. 


Equation 
—-— Equation (62) 
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Pressure coefficients C, vs. Mach number cal- 
culated from Eqs. (62) and (65). 
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However, with increasing Mach number the local ve- 
locity at the surface increases faster than in the undis- 
turbed stream and this causes an increase of the skin 
friction. Another effect is due to the increase of ad- 
verse pressure gradient at the rear portion of the air- 
foil causing an increase in thickness of the boundary 
layer. The increase of the adverse pressure gradient is 
caused by the large increase of the maximum suction 
shown in the sections on Method of Small Perturbations 
and Applications of the Hodograph Method. Both ef- 
fects can be taken into account—at least approximately 
—by the theoretical considerations given in the sec- 
tion on Method of Small Perturbations. It was pointed 
out that the velocity and pressure distribution along the 
surface of a slender body or an airfoil at a certain Mach 
number is approximately the same as in a flow of an in- 
compressible fluid provided the body is expanded later- 
ally in the ratio 1:+/1 — M,°. This reasoning furnishes 
a method for estimating the compressibility effect on 
the drag in the following way. Consider first, sym- 
metrical sections. If the thickness ratio of the airfoil 
is denoted by ¢, the fictitious thickness corresponding to 
the Mach number M, is equal to t/1/1 — M2. Pro- 
vided the influence of the thickness on the drag is 
known for the family of airfoils considered, the mini- 
mum drag coefficient Cp, can be put approximately 
equal to the minimum drag coefficient Cp, of a section 


of the thickness ¢/+/1 — M,*. To estimate the correc- 
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tion for an arbitrary angle of attack it has to be taken 
into account that by the lateral expansion of the flow 
pattern the angle of attack is also increased in the ratio 
1/\/1 — M,2. Therefore, the drag coefficient of a sec- 
tion of thickness ¢ at the angle of attack a and Mach 
number MM, will be approximately equal to the drag 
coefficient of a section of the thickness t/+/1 — M,? at 
the angle of attack a//1— M2. If the effect of com- 
pressibility on the lift coefficient is known, one can use 
the drag coefficient of the expanded section at the lift 
coefficient Cz ,,, which corresponds to the Mach num- 
ber \/, at the angle of attack a. For cambered sec- 
tions the camber has to be also increased in the ratio 


1:1 — M2. 
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N.A.C.A. 4412 at a = 1°52’. 

Fig. 14 gives a comparison between values of the drag 
coefficient of the section N.A.C.A. 4412 measured at 
1°52’ angle of attack and computed according to the 
method explained above. The influence of thickness 
was evaluated according to N.A.C.A. Technical Report 
No. 669 which apparently represents the last word in 
this somewhat controversial question. It was assumed 
that the relative influence of the thickness is independ- 
ent of Reynolds number. It is seen that the accordance 
between the measured and estimated values is very 
good, until the “compressibility burble’’ occurs, indi- 
cated by simultaneous sudden drop of the lift coefficient 
and increase of the drag coefficient. , 

The ‘‘compressibility burble’’ is apparently caused by 
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a shock wave extending from the surface to a certain 
distance from the surface. Since the magnitude of the 
pressure or the velocity jump is variable along the dis- 
continuity surface, the flow behind the shock wave is a 
flow with continuously distributed vorticity even if the 
flow is irrotational in front of the discontinuity surface. 
Since far behind the airfoil as the flow becomes parallel 
again—at least theoretically—the shock wave must 
cause a non-uniform velocity distribution; 7.e., a wake 
in addition to the wake caused by the boundary layer. 
This has been experimentally shown by J. Stack, W. F. 
Lindsey, R. E. Littell and also by C. N. H. Lock. 
There are only a few measurements available for drag 
and lift of airfoils beyond the velocity of sound and it 
seems that nobody has been able yet to obtain reliable 
measurements at the speed equal to the velocity of 
sound. The observations of the N.A.C.A. are extended 
essentially to the speed at which the compressibility 
burble occurs. The experimenters working with the 
high-speed tunnel in Guidonia found by increasing the 
speed beyond this limit a ‘‘zone of oscillation” in which 
no stable equilibrium can be established between the 
air forces and the balance. Beyond this zone the drag 
coefficient has a higher value and is apparently con- 
stant in the range of the velocity of sound. Fig. 15 
shows the Italian measurements for an airfoil with 6.32 
per cent thickness whose shape is almost identical with 
that of a Clark Y section. Curiously enough the 
N.A.C.A. measurements on a Clark Y airfoil of 6 per 
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Drag coefficient Cp vs. Mach number M, at 
Cr, = 0.45. 

(I) Italian propeller section of 6.32 per cent thick- 

ness (A. Ferri). 


Fic. 15. 


(II) Clark Y section of 6 per cent thickness (J. 
Stack). 
(III) Géttingen 622 section of 8.85 per cent thickness 


for a = 2° and 4° at M, = 1.47 (A. Buse- 
mann and O. Walchner). 
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Fic. 16. Drag coefficient Cp of R.A.F. 6 section, 10 


per cent thickness, at a = 2°. 
(I) N.A.C.A. tests (J. Stack). 
(II) British propeller tests (G. P. Douglas). 
(III) Open jet tests (L. J. Briggs and H. L. Dryden). 


cent thickness are at wide variance with the Guidonia 
results. The two curves correspond to equal values of 
C, ~ 045 for M, — 0. For comparison the drag 
coefficients of a similar airfoil of 8.85 per cent thickness 
at M, = 1.47 for 2° and 4° angle of attack are indicated. 
These data are taken from measurements of A. Buse- 
mann. The N.A.C.A. data and results of British pro- 
peller tests concerning the compressibility effect on the 
drag of a 10 per cent thick R.A.F. 6 section at 2° angle 
of attack are in fair accordance (Fig. 16). In the same 
diagram the data obtained by Briggs and Dryden in an 
open high-speed jet are also plotted. It appears that 
these values are probably too low. 

There are some measurements on lift and drag of a 
section with circular-are upper and flat lower surface 
between MM, = 1.47 and 2.13. These data are plotted in 
Figs. 17 and 18. In order to complete the picture as 
far as possible the corresponding data for the N.A.C.A. 
4409-34 section are plotted for 0.35 < M, < 0.84. 
This section is somewhat similar, having its maximum 
thickness at 40 per cent of the chord (instead of 50 per 
cent) and an almost flat bottom. It seems that the 
general shape of the drag curve is similar to that of the 
drag curve of a sharp-nosed projectile. The lift ap- 
parently decreases greatly after the first occurrence of 
the compressibility burble. The rise shown by the last 
points needs further confirmation. The theory of lift 
and drag of airfoils for supersonic speed has been de- 
veloped by J. Ackeret, A. Busemann and S. G. Hooker. 
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Fic. 17. Lift coefficient Cy vs. Mach number M. 
(I) N.A.C.A. 4409-34 (J. Stack and A. E. v. 
Doenhoff). 
(II) Circular-are section of 8.85 per cent thickness 
(A. Busemann and O. Walchner). 
(III) Circular-are section of 8.8 per cent thickness 
(A. Ferri). 


There are no data published on the compressibility 
effect on the drag of streamlined bodies such as fuse- 
lages and nacelles. In general it can be expected that if 
a certain streamlined shape has a favorable drag coef- 
ficient at low speed, it has to be expanded in the length 
direction or contracted laterally in a linear ratio about 
equal to 1:\/1 — M,? in order to keep the drag coef- 
ficient low at high speed. 


THE CRITICAL MacH NUMBER 


The Mach number at which the compressibility 
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Fic. 18. Drag coefficient Cp vs. Mach number M. 
(I) N.A.C.A. 4409-34 (J. Stack and A. E. v. 
Doenhoff). : 
(II) Circular-are section of 8.85 per cent thickness 
(A. Busemann and O. Walchner). 
(III) Circular-are section of 8.8 per cent thickness 
(A. Ferri). 
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burble first occurs is called the critical Mach number. 
In general it is defined as the Mach number of the un- 
disturbed flow for which the local velocity at some point 
of the surface reaches the local velocity of sound. Us- 
ing Glauert’s approximation the condition for the criti- 
cal Mach number is given by 


V+,’ +a (66) 


Assume that the pressure coefficient of a given airfoil or 
slender body at the point of maximum suction is equal 
to —C,; then the corresponding pressure coefficient 


for the Mach number M, is 
“Cy. = —C,/V1 — M? 
and since 
Cy, = 20,’/V 
the Eq. (66) can be written in the form 
V1 + (Cp, /2)] = a1 


or 


mf + ( C;,/2V1 = M:) | mo (67) 
Hence, one obtains the following equation for the criti- 
cal Mach number 

2[(1 = M;)'“(1 + M,)'”] /M, = 


Cy, (68) 


E. Jacobs, by using Eq. (8) for the relation between 
pressure and velocity, obtained the expression: 


C,, = 


——. 2+ — 1)MA7/0-Y 
oV1 — Me E a ( oss ae | [ae 


7? 1 
(69) 





However, the question arises whether it is justified to 
introduce expressions involving higher terms in a linear- 
ized theory that neglects squares and products of the 
perturbation velocities. Possibly the error made in 
evaluation of the ratio C a! Cp, is of the same order as 
the correction involved by use of the exact thermody- 
namic relation (8) instead of (66). 

The method developed by Tsien and the author per- 
mits the calculation of the critical Mach number by 
equating the pressure coefficient computed by Eq. (62) 
with the value of the pressure coefficient corresponding 
to a local Mach number M = 1 according to Eq. (8). 
In Figs. 10 to 12 this limit is plotted by a full line. 

Fig. 19 shows a comparison of the relation ‘‘critical 
Mach number’’ vs. C, according to the different ap- 
proximate theories. For comparison a few experimen- 
tal results are also indicated. 

Of course, by laborious analytical methods it would 
be possible to refine the calculation of the critical Mach 
number, defined by the first coincidence hetween local 
velocity and local velocity of sound. However, it ap- 
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pears questionable whether this is the right criterion 
for the first appearance of a compressibility burble. A 
closer examination of the fundamental mathematical 
problem shows that at least in an ideal compressible 
fluid a shock wave does not necessarily occur if the 
local velocity exceeds the velocity of sound. The hodo- 
graph method can be used to throw some light on this 
question. 

It was shown in the section on the Hodograph 
Method, that if ¥(v, 8) is the stream function for a flow 
pattern of a compressible fluid where v is the magnitude 
of the velocity and 0 the angle of inclination of the veloc- 
ity vector at an arbitrary point, then y as a function of 
v and 0 must satisfy the differential equation 


1 — M? dy (° oy’ 
EE. Ba ages = 0 
p> oe? p Ov v 


Now it is possible to find exact particular integrals 
of this differential equation, and such solutions are sup- 
posed to give exact flow patterns of a compressible ideal 
fluid. The physicist or engineer is, of course, not satis- 
fied by having determined streamlines in the hodograph 
plane and is anxious to see the flow pattern in the physi- 
cal x, y plane. If the function ¥(v, 6) is known, the 
computation of streamlines ¥(x, y) = const. in the 
physical plane requires a further integration. How- 
ever, in performing this integration it is found that in 
some cases a streamline that is perfectly continuous in 
the hodograph plane shows a singular behavior: plot- 
ting successive points of a streamline in the physical 
plane at a certain point which can be called the critical 
point, the streamline cannot be continued in the flow 
direction, and gives a quite impossible flow pattern with 
two families of crossing streamlines in the same plane. 
The velocity at this critical point is neither zero nor in- 
finite. However, by continuing the calculation one ob- 
tains a second branch of the streamline, which meets 
with the first branch in a kind of ‘‘cusp.”’ It appears 
that the flow cannot enter in a certain “forbidden” 
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Fic. 19. Critical Mach number M, vs. maximum suction- 
pressure coefficient at low speed C),,. 
(I) Eq. (68). 
(II) Eq. (69) (E. Jacobs). 
(III) Hodograph method (Th. von K4rm4n and 
H. S. Tsien). 
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region without violating one of the physical assumptions 
underlying the analysis. The two assumptions are: 
continuity and irrotationality. Since the continuity 
cannot be violated, it must be assumed that the flow be- 
comes rotational. In an ideal compressible fluid this 
can happen only through a pressure jump; 7.e., a shock 
wave. In this sense the occurrence of such critical 
points indicates the necessary appearance of shock 
waves. 

The phenomenon of critical points was found first 
by M. and F. Clauser who, on the suggestion of the 
present author, calculated a family of exact solutions of 
the Eq. (70) and discussed the corresponding flow pat- 
terns. They investigated especially the flow bounded 
by two inclined straight walls. Recently, F. Ringleb, 
by a method similar to that employed by the Clausers, 
found exact flow patterns of compressible fluids. He 
presents in his paper one example which is especially in- 
structive as to the question of the first appearance of 
shock waves and will be briefly discussed here. 





° 


Fic. 20. Flow in a curved convergent-divergent nozzle 
(F. Ringleb). 


Fig. 20 shows a kind of nozzle which is essentially a 
curved convergent-divergent nozzle. The boundaries 
are exact streamlines of a compressible fluid obtained 
in the manner explained above. The velocity is zero at 
the entrance and at the exit at infinity. The lines of 
constant velocity—which happen to be exact circles— 
show that the flow is partly subsonic, partly supersonic. 
Consider first the nozzle bounded by the lines (a) and 
(6). The velocity distribution along the two walls is 


353 


given by the lines (a) and (0) in the lower part of Fig. 20. 
It is seen that the flow is perfectly smooth and continu- 
ous. However, if the nozzle is bounded by the lines 
(a) and (c) the velocity has infinite slopes at the points 
P and Q. If one plots another streamline proceeding 
forward inside, it is found that the streamline stops at a 
certain point and turns back, crossing the nozzle; this 
is, of course, physically impossible. 

It is remarkable that the maximum velocity reached 
at the wall (c); 7.e., in the limiting case, is as high as 
1.67a9 (where dp is the velocity of sound at rest) and that 
the velocity at the point where a discontinuity first ap- 
pears in the flow is equal to 1.20ap. 

The German author does not discuss further the con- 
ditions for the singular behavior of the streamlines. 
However, the equations developed in the section on the 
Hodograph Method of this paper make it possible to 
determine a condition for the occurrence of an infinite 
acceleration. Obviously the acceleration in the direc- 
tion of the flow is equal to (0v/0s)(ds/dt) = (Qv/ds)v = 


(Ov/O¢g)v?. Hence, if v has a finite value, the accelera- 
tion can be infinite only if0v/Og = ©. Nowaccording 
to Eq. (48) 


Ov/Og = —(1/A)(Oy/00) 


where A = (0¢/0v)(OY/00) — (0¢/00)(Oy/dv) is the 
so-called Jacobian previously discussed. If the case 
Oy /00 = o is excluded for the time being, dv/Oy can be 
infinite only when 


A = (0¢/dv)(OY/00) — (Og/00)(OY/Ov) = 0 (71) 


Substituting the values of 0g/dv and 0¢/00 from the 
Eq. (50) one obtains the condition 


(1 — M*)y/00)? + v°(Oy/dv)? = 0 


The expression on the left side is obviously always posi- 
tive if M < 1. Therefore, the condition (72) can be 
only satisfied if M > 1; 1.e., if the flow is supersonic. 
Eq. (72) has a simple geometrical meaning. Denote 
the angle of inclination between a circle v = const. and 
the tangent to the line y = const. at an arbitrary point 
of the hodograph plane by 8. Then 


(72) 


1 Oy\ /dy ‘ 
=—(-—]) /— 3 
tan ( 2) / (73) 
Hence, from Eq. (72) follows 
tan? 8B = 1/(1 — M?) (74) 


In the theory of supersonic flow the angle sin—! (1/M) 
= ais called Mach’s angle; it is equal to the angle be- 
tween the velocity vector and the so-called stationary 
wave-front produced by a small perturbation applied 
to the fluid. Hence, according to Eq. (74), 8 = ta. 
The lines intersecting the circles v = const. at the angle 
a are called the real characteristics of the differential 
Eq. (70) and for given constant y are epicycloids. 
Hence, infinite acceleration can occur at points where 
the curve yy = const. in the hodograph plane is tangen- 
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Fic. 21. Condition of infinite acceleration and breakdown 
of isentropic irrotational flow. 


tial to a characteristic curve of the differential Eq. 
(70)* (Fig. 21). 

The direction of the tangent to a y = const. curve 
indicates the direction of the infinitesimal change of 
the velocity vector; 7.e., the direction of the resultant 
acceleration, in the sense that the angle between the 
velocity and acceleration vectors is the same in the 
hodograph and physical planes. Hence, at the point 
where infinite acceleration occurs the angle between the 
resultant acceleration and the normal to the velocity 
is equal to the Mach angle. Applying this result to 
the case that the velocity along a streamline changes 
from subsonic to supersonic and back to subsonic again 
one can show that in order to avoid infinite acceleration 
the angle between the acceleration and the normal must 
be smaller than the local Mach angle in the entire 
supersonic range. Further investigations are necessary 
to expand the physical meaning and importance, if 
any, of this finding. 

In Fig. 22 streamlines for the flow pattern shown in 
Fig. 20 are plotted in the hodograph plane. Fig. 23 
refers to the velocity distribution found experimentally 
for the N.A.C.A. 4412 section at various Mach num- 
bers. The lines are plots of the “zero streamline” 


* The fact that the singular behavior of the streamlines occurs 
at points where streamline and characteristics in the hodograph 
plane have the same direction was first found by the Clausers. 
The condition (71) is mentioned in the paper of Ringleb. It ap- 
pears that if a streamline is tangent to a characteristic in the 
hodograph plane, in the physical plane the characteristics (wave- 
lets) of the same family meet; their envelop probably introduces 
the first discontinuous wave front. 
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Fic. 22. Streamlines and characteristics in hodograph 
plane for a curved convergent-divergent nozzle. 


(upper boundary) in the hodograph plane; in other 
words the velocity vectors occurring along the boundary 
are plotted as to magnitude and direction and their end 
points connected by a curve. It is seen that the line 
at a point close to the point where the shock wave oc- 
curs is tangential to characteristics passing through the 
same point. However, this coincidence might be only 
accidental because the velocity distribution is greatly 
changed by the shock wave, whereas the theorem proved 
above refers to the first occurrence of an infinite slope 
of the pressure, but not to a finite pressure jump. 

Flow patterns with continuous transition of the 
velocity of sound in both directions have been calcu- 
lated before by G. I. Taylor, W. Tollmien and recently 
by H. Gértler who investigated the flow along a wave- 
shaped solid wall by computing a second approximation 
on the lines of the Glauert-Prandtl method. 

The author does not wish to exaggerate the impor- 
tance of the above considerations. It is merely shown 
that in the case of an ideal compressible fluid continu- 
ous, shockless flow patterns can be constructed in which 
the velocity of flow reaches magnitudes as great as 1.6 
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times the velocity of sound and becomes subsonic again. 
However, this does not prove yet that such flow pat- 
terns are really possible in a real fluid. Furthermore, 
it is possible that the flow with shock wave is more 
stable than the continuous flow pattern. If a body 
moves with supersonic velocity in a compressible fluid, 
shock waves appear unavoidable, 7.e., flying with super- 
sonic velocity will always involve wave resistance in 
addition to the frictional drag. However, it appears 
that careful theoretical and experimental research 
might be able to push the velocity of flying nearer to 
the velocity of sound than is possible now. The mere 
fact that the air passes over a wing with supersonic 
velocity does not necessarily involve the occurrence of 
a compressibility burble and energy loss by shock wave. 
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Letter to the Editor 


June 12, 1941 


Dear Sir: 

The following discussion applies to the article, Permanent 
Buckling Stress of Thin-Sheet Panels under Compression, by 
Messrs. W. L. Howland and P. E. Sandorff published in the 
May, 1941, issue of the Journal. 

On page 264, the authors derive a value for the strain at the 
defined yield strength (Equation (11)) and subsequently use this 
strain as the maximum allowable total strain for the panel. The 
value which the authors use is apparently based on tensile-test 
data for 24ST Alclad sheet. Long before this strain is reached, 
the Alclad coating begins to yield. (Yield strength of com- 
mercially pure aluminum is 4000 Ibs. per sq.in.; see Structural 
Aluminum, Aluminum Company of America, 1938, page 10.) 
In a tensile (or compression) test, the core material, because 
of its greater thickness, has a predominating influence. There- 
fore, if the core material has not yielded, the residual strain in 
the Alclad sheet will be small. When initial bending of the 
sheet is involved, as in buckling, the yielding of the Alclad coat- 
ing, which is a large distance from the neutral surface, is more 
important than in the case of the tensile-test, and when per- 
manent bending is involved, as in the present case, the yielding 
of the Alclad coating has a tremendous effect. (This has recently 
been checked by tests at the research laboratories of the Alumi- 
num Company of America.) It is therefore believed that the 
value of maximum strain based on tensile test data has little 
significance in connection with permanent buckling of Alclad 
sheet, which depends wholly on the bending characteristics of 
the sheet. Since the authors subsequently (page 265) chose an 
empirical value for the maximum strain, the preceding comments 
on the action of Alclad sheet are not too important in connection 





with the value used for the maximum strain. However, the fact 
that the Alclad covering yields at low strains is of utmost im- 
portance in interpreting the buckle-depth data. 

Referring to Fig. 5 on page 265, when the specimen was first 
stressed to 5500 Ibs. per sq.in., the total strain was great enough 
to cause yielding in the Alclad coating. When the load was 
released, the bending strength of the unyielded core was suffi- 
cient to return the sheet practically to its original shape despite 
the slightly yielded Alclad coating. However, when the critical 
stress was approached in the next loading, the effect of the pre- 
stressed Alclad coating, being magnified by the presence of the 
compressive stress, appreciably increased the buckle depth. 
Only when the stiffener stress reached the previous maximum 
of 5500 Ibs. per sq.in. was the buckle depth the same as the 
buckle depth during the first loading. After continuing the 
loading up to 7800 Ibs. per sq.in. and then releasing, the core 
again returned the sheet practically to its original shape. On 
the next loading, the Alclad coating, having more residual stress 
than on the previous loading, caused still deeper buckles in the 
vicinity of the critical stress (below previous stiffener stress). 
On subsequent loadings, to higher stresses, the increased yielding 
of the Alclad coating was sufficient to cause measurable per- 
manent buckles at the first increment of load. With this ex- 
planation of the data plotted on Fig. 5, the authors’ conclusion 
(see page 268), ‘‘That an ideal (absolutely flat) panel once sub- 
jected to stress greater than the critical-buckling stress will 
exhibit some residual roughness at any lower reference stress, 
even at no load,” is not ‘difficult to believe.’’ 

Wiey P. MonTGOMERY 
Vought-Sikorsky Aircraft Division, 
United Aircraft Corporation 
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Fog Study of Miami, Florida 


JOHN RANDLE DEeHART 


Eastern Air Lines, Inc. 


As advection plays an important role in every case of 
low visibility for this area, the geographical location of 
Miami is of great importance. Miami has an extreme 
water exposure, located as it is, with the Atlantic Ocean 
on the east and with the Gulf of Mexico approximately 
75 miles to the west. To the northwest is Lake Okee- 
chobee and the adjacent Everglades. 

The Everglades contribute not only moisture to the 
lower layers of air, but also at frequent intervals add 
quantities of smoke from burning muck. This may or 
may not be a contributing factor to fog formation, but 
certainly it does occasionally cause visibilities to become 
so low as to be dangerous to aircraft. In fact, one-third 
of Miami’s visibilities of one mile or less, at surface 
winds of less than 5 miles per hour, are due to smoke. 

In view of the topography, almost any air arriving at 
Miami, regardless of its path, must have had some re- 
cent water history. The graphic story of this is told in 
Fig. 1, where it can be seen that air, arriving over trajec- 
tory A, has its path almost entirely over the Everglades, 
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Fic. 1. Topography and trajectories. 


while air arriving over trajectory B has had some history 
over the Gulf of Mexico, and only the latter part of its 
journey was over the Everglades. Air in both trajec- 
tories A and B will become stabilized during the night 
when the air temperature is higher than the temperature 
of the water of the swamp land. 


Awarded first prize for 1940 by the Air Transport Association 
of America “for the best research, by an airline employee, having 
practicable application to meteorology or dispatching.” 

This is the first such award given by the Air Transport Associa- 
tion of America. 
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Fic. 2. Distribution of fog by gradient direction. 

Air in trajectory C, which has a high moisture content 
due to its water history over the Caribbean, is stabilized 
by its passage over the progressively cooler shore waters. 
Later when it passes over the city, it is exposed to 
smoke and dust pollution. These two factors tend to 
cause the largest amount of fog of any other gradient 
direction. 

In consideration of trajectory D: This air is inher- 
ently unstable because of its passage over the warm 
waters of the Gulf Stream, the edge of which is approxi- 
mately 12 miles from the airport. In this short path 
the air has little opportunity to cool and become stable, 
except with the very lightest of winds. 

A gradient windrose shows air arriving from NE to E 
produces the least amount of fog of any other gradient 
direction, and that from E to SW produces the greatest 
amount; the first because of this unstabilizing effect of 
the Gulf Stream, and the second because its path has 
been parallel to the sea isotherms of the Gulf Stream 
and during the latter part of its history it has become 
stabilized over the cooler bay waters. (See Fig. 2.) 
Fog occurs with winds from SW to W more often than 
from W to NE because this air is usually more moist due 
to its water history over the Gulf of Mexico. 

Compare this gradient windrose with the surface 
windrose of restricted visibility as given in Fig. 3a, 
which shows a maximum of low visibility from NW. 
However, this surface chart includes in it low visibility 
due to smoke and rain, while the gradient windrose 
excludes both of these. The light gradients from SE to 
S would most likely show up as calms on the surface 
windrose. 
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Distribution of fog by surface winds and distribution of 
fog by months. 


Fic. 3. 


SEASONAL VARIATIONS 


Another consideration in this study is the seasonal 
variation of low visibility. (See Fig. 3b.) The maxi- 
mum amount of fog in this area occurs during the month 
of March, which only slightly exceeds the fog formation 
during January. The drop in the fog curve for Febru- 
ary may be explained in that this month has the highest 
average wind velocity, and the fact that low visibilities 
occur less frequently from May through September is 
due to the greater amount of insolation at this time. 

In breaking the cases of fog down into their respec- 
tive air masses the following result was found, which is 
about the normal occurrence of these air masses: 


43.5 per cent 
38.7 per cent 
17.8 per cent 


Tropical Maritime 
Modified Polar Continental 
Modified Polar Pacific 


Therefore, it can be seen that fog may form in any 
air mass that covers Miami. Stratus occurs very 
rarely at Miami, if at all, unless in conjunction with a 
front. These frontal cases have been discussed fully 
in such papers as George’s The Causes and Forecasting 
of Low Ceilings and Fogs at Jacksonville Airport (pub- 
lished by Eastern Air Lines, Inc., 1939), and will not be 
restated here. 


QUALITATIVE ESTIMATION 


Whether or not fog will form, and if so, will it be 
light or dense, may be computed from the curves 
shown in Fig. 4, which are given for both Tropical and 
Polar Air Masses. These curves are not to be used if 


local showers have occurred in the four hours preceding 
the 1930 E map because the dew point would then be 
fallaciously high. To use these curves, plot gradient 
wind as given by isobars on the 1930 E map against the 
dew point depression as given on that same map. If 
the point so plotted falls between the inner and outer 
curves, fog will occur, but the visibility will not be be- 
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low one-half mile; but if it lies within the inner circle 
the fog which occurs will cause the visibility to be be- 
tween zero and one-half mile. 

The gradient wind must also be forecast for the early 
morning hours because a point that plotted in the light 
fog sector might move into the dense fog sector if the 
gradient wind should decrease; the reverse might be 
true if the gradient wind should increase. 


TIME ESTIMATION 


After deciding that fog will form, the next step is to 
determine the time of fog formation. Using the curve 
of Fig. 5, plot the sunset dew point depression on the 
curve and read the time of formation from the abscissa. 
The curve assumes constant gradient winds between 9 
to 16 m.p.h. For decreasing winds fog forms 1 to 3 
hours earlier tian indicated; and for increasing winds, 
1 to 3 hours later. This curve, of course, is to be used 
when the forecaster expects the cloud cover to dissipate 
soon after sunset. If this is not the case, suitable cor- 
rections must be made. 

To correct for the presence of cloud cover, the fore- 
caster must add 75 per cent of the number of hours 
that the cloud cover is expected to remain, to the value 
derived from the time of formation curve. 

The curves themselves are empirical and were de- 
rived from plotting U.S. Weather Bureau data from 
Eastern Air Lines weather charts. 


FORECASTING RADIATION-ADVECTION Foc 


Presentation of the following two cases, (1) a normal 
radiation-advection fog, and (2) a case when fog did 
not form, show how the forecast was simplified and 
more accurate results were obtained by the use of 
charts and curves compiled for this study. 

Ground fog may be expected in Miami in the normal 
case shown in Fig. 6, in view of the facts: 


1. Isobars show gradient wind from 110° at 15 


m.p.h. 
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Fic. 4. Qualitative forecasting curves. 
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Fic. 6. Normal case. 














Fic. 7. Case when fog did not occur. 


























Fic. 9. Passing of front. 
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Fic. 10. Lifting and breaking of fog. 
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2. Decreasing gradient winds are expected because 
of rising tendencies at Key West and Fort 
Meyers. 

3. Gradient of moisture shows no advection of drier 
air. 

4. Dew point at sunset in Miami was 64°F., and the 
temperature 74°F. 


With a dew point as high as 64°F., Miami is un- 
doubtedly in Tropical Maritime Air. Plotting gradient 
wind against dew point depression at 1930 E on Fig. 4 
for Tropical Maritime Air, the point falls in the light 
area. But, because the gradient wind is decreasing, the 
point will move from the light to the dense fog area. 

The result of plotting the temperature and dew point 
depression on the curve of Fig. 5 shows that fog should 
form at 0200 E. Actually, at 0135 E visibility was re- 
duced to one-half mile. By 0235 E visibility was reduced 
to one-fifth mile in dense ground fog, and remained 
between one-eighth and one-half mile until 0700 E. 

The importance of forecasting and of noting all de- 
tails such as gradient of moisture, forecasting gradient 
wind and trajectory cannot be overlooked, because the 
history of the air arriving at Miami can be derived only 
from the few stations that are located on the Florida 
peninsula. As, for example, in Fig. 7, using the gradi- 
ent wind from the map, with a temperature of 67°F. 
and a dew point of 60°F., without further investigation 
it would appear that dense ground fog would form. 
However, employing the graphical method of forecast- 
ing the gradient wind, using an isallobaric component 
as shown on Fig. 7, it is clear that the gradient wind 
should be 160° at 11 m.p.h. by 0130 E, and 190° at 
19 m.p.h. by 0730 E. From this increasing gradient it 
is apparent that any fog forming would be light. Fur- 
ther study shows dew points of 52°F. at Key West and 
54°F. at Fort Meyers; therefore, any advection of air at 
Miami would be drier air, further reducing the possi- 
bility of dense fog. 

Actually, the balloon runs gave these values: 


1800 E 110° 5m.p.h. 
2500 E 180° 10 m.p.h. 
0600 E 230° 16 m.p.h. 


and the visibility did not go below two miles at any 
time during the night. The following morning showed 
that the dew point had dropped to 51°F. Such accu- 
rate computations of wind acceleration are, of course, 
possible only with unusually smooth isallobars. 


FORECASTING Post COLD FRONTAL FoG 


Another type of fog peculiar to this area is the post 
cold frontal fog brought about by two tropical air 
masses. The following sequence of events, of which 
this is a typical example, produces the worst cases of fog 
and low ceilings at Miami. (See Fig. 8.) The Tg air, 
shown as trajectory A, has had an approximate three- 
day history over water. It moves eastward in the 


circulation brought about by the advancing cold front 
shown in the upper left-hand corner of the map. The 
Tg air underruns the very warm Ta air, and a NE-SW 
line of frontogenesis is formed. This frontogenetic line 
continued southward and crossed Miami near 2300 E. 

Fig. 9 shows the situation 6 hours later. The line of 
frontogenesis has now become a true 7a-7g front, and 
the original cold front is undergoing frontolysis. At 
Miami the stabilizing effect of the warm air aloft, plus 
the rapid clearing behind the front coupled with low 
wind velocities, and perhaps some mixing, caused visi- 
bility to be reduced to one-half mile by 0335 E, and to 
one-eighth mile by 0635 E. This fog lifted to stratus 
due to the advection of unstable air brought about by a 
wave formation off the coast as shown in Fig. 10. This 
same condition lowered ceilings to zero and visibilities 
to one-eighth mile at Key West, a most unusual occur- 
rence. 

The author feels that careful consideration and 
weighing of all available data should result in the rea- 
sonably accurate fog forecasting in at least 80 per cent 
of fog cases, and that errors in forecasting are directly 
attributable to the lack of sufficient time for proper 
consideration, to inadequate data, or to carelessness on 
the part of the forecaster. 


CONCLUSIONS 


Combination of conditions under which fog will 
form and visibility will lower to one-half mile or below: 


1. Forecast from 1930 E map. 
2. Showers have occurred two to four hours pre- 


vious. 

3. Gradient wind will not accelerate to 20 m.p.h. by 
0730 E. 

4. Dew point depression on 1930 E map is not more 
than 11°F. 

or 
1. Forecast from 1930 E map. 
2. Showers have not occurred two to four hours 


previous. 

3. Gradient wind is below 10 m.p.h. and shows de- 
celeration, or at least no acceleration. 

4. Dew point depression on 1930 E map is not more 
than 5°F. 


Combination of conditions under which fog will form 
and visibility will lower to one-half to two miles: 


1. Forecast from 1930 E map. 

2. Dew point depression is between 5°F. and 11°F. 

3. Gradient wind not less than 10 and not more 
than 20 m.p.h. up to 0730 E. 


This study was made from Eastern Air Lines’ maps, 
drawn and analyzed at six-hour intervals from 1935 to 
1940. All frontal analyses and air mass designations 
are their own. 
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The Use of Liquid Oxygen for High- 
Altitude Flying 


JOHN D. AKERMAN 


University of Minnesota 


SUMMARY 


A method and apparatus is described to produce gaseous oxygen 
for use in high-altitude flying, from liquid oxygen. The elimina- 
tion of the danger of discontinuity of oxygen boiling and provision 
for securing the required volume of oxygen gas is accomplished by 
keeping the liquid oxygen under a predetermined pressure and 
supplementing the normal evaporation by spilling a necessary 
amount of liquid oxygen in the gas contained. Fast evaporation 
takes place to maintain a constant supply of oxygen gas inside the 
apparatus which is connected with a second-stage reducing valve 
of a B-L-B oxygen mask as developed by the Mayo Foundation. 

To secure the spilling of liquid oxygen into the gas container, 
an absolute pressure valve is provided in the gas contained to 
keep a predetermined pressure of gas in the liquid container. 
The liquid container does not require the use of a vacuum bottle. 
Although no attempt is made to design the apparatus for any 
specific airplane, the mechanism is flexible for adaptation to meet 
any specific requirements of oxygen supply. Apparatus for han- 
dling and storing liquid oxygen and the possibilities of a portable 
liquid oxygen producing plant are discussed. 


INTRODUCTION 


a OXYGEN as a source of supply of oxygen gas 
for high-altitude flying was used extensively by the 
Military Service from 1925 to 1934. Since that time, 
for convenience and safety, the use of compressed oxy- 
gen gas has been introduced. Why the use of liquid 
oxygen was abandoned is not a matter for discussion 
in this paper, but there is sufficient evidence that with 
improper liquid oxygen installations there is a danger 
of shortage of oxygen gas during flight when descending 
from higher to lower levels. This fact was brought out 
by Dr. Jean Piccard* whose conclusions were: 

During ascent the boiling point of oxygen goes down 
due to decrease of pressure. On descent, due to in- 
crease of pressure, oxygen stops boiling. At this period 
the pilot will be short of oxygen. There are several 
cases on record, here and abroad, where such a situation 
was encountered by pilots using liquid oxygen installa- 
tions. 

After discussing with the author the possibilities of 
forced vaporization of oxygen under controlled pressure, 


Presented at the Physiological Problems session, Ninth Annual 
Meeting, I.Ae.S., New York, January 29, 1941. 

Charles Cole, Research Assistant, University of Minnesota 
Engineering Experiment Station, Minneapolis, Minnesota, and 
Norvin Erickson, Research Assistant, Mayo Foundation, 
Rochester, Minnesota, aided in the construction of apparatus 
and in the collection of data. 

* Change of Boiling Point of Liquid Oxygen at High Altitudes, 
Dr. Jean Piccard. 


Dr. Piccard fully endorsed this system of using liquid 
oxygen as supply for high-altitude flying. No doubt, 
the problem would not have reached its present stage 
without his help. 
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Fic. 1. Small size oxygen gas 
cylinder capacity (1 liter of liquid 
oxygen = approximately 880 liters 
of gaseous oxygen). 


Fig. 1 shows a compressed oxygen gas cylinder suffi- 
cient for one man for 1 hour and 59 minutes at 30,000 
ft. altitude if permissible pressure is 1800 lbs., but if the 
pressure is reduced to 400 Ibs. it is sufficient for only 26 
minutes at 30,000 ft. The decrease of pressure from 
1800 Ibs. to 400 Ibs. decreases the contents from 203 
liters to 45 liters. One liter of liquid oxygen would 
produce approximately 880 liters of gaseous oxygen at 
sea level and would be sufficient for one man for ap- 
proximately 8'/2 hrs. at 30,000 ft. Installation with 
1800 Ibs. pressure was satisfactory until it was dis- 
covered that if such a cylinder is hit by a bullet, an 
explosion occurs which is very damaging to the air- 
craft. 

Therefore, specifications are now in force to make 
new installations of compressed oxygen on military air- 
craft in containers with a maximum pressure as low as 
400 Ibs. per sq.in. Such a container for the same 
amount of oxygen as shown in Fig. 1 would be 5.1 times 
larger and presents considerable problems of installation 
on modern aircraft. There is also the problem of pro- 
viding armor against machine gun bullets. 

To introduce again the use of liquid oxygen in order 
to save weight and space, and to meet requirements for 
large volume of oxygen, seems very promising except 
for two major difficulties. The first and easiest is 
technical development, and the second, and hardest, 
prejudice created by past experience. 

This paper will be devoted to the first item, leaving 
the combat of prejudice to the merits of the develop- 
ment. Before going into details of the method and 
equipment, it is essential to make clear the basic condi- 
tions governing this work. 
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1. No attempt is made to produce an actual instal- 
lation for any airplane, commercial or military, 
because requirements and conditions vary so 
that it is impossible for educational institutions 
with limited financial means, to meet any con- 
ditions specifically. Specific adaptation and 
expensive refinements can be left for the 
future. 

2. Requirements for the amount of oxygen for one 
man at different altitudes used in this work are 
those supplied by the Aviation Medicine Divi- 
sion of the Mayo Foundation. 

3. Equipment is designed to be usable directly with 
the B-L-B mask developed by the Mayo 
Foundation. 

4 Accepted as standards are flow meters as used by 
the Mayo Foundation since some other equip- 
ment has shown dangerous variations at high 
altitudes. 

5. Difficulties of producing, handling and storing 
liquid oxygen may be in the future overcome 
with technical developments described later 
in this paper. 


The first requirement in the use of liquid oxygen is 
the elimination of the possibility that there will not be 
enough oxygen gas delivered at some time during the 
flight. 

The second requirement is the elimination of the 
possibility of building up high pressure inside of the 
container which might then burst and release oxygen 
gas at high pressure. 

There is an unfounded fear that spilled liquid oxygen 
or oxygen gas at low pressure is a dangerous explosive 
by itself. 

Fig. 2 shows a diagram of the setup in connection 
with the B-L-B mask with standard second-stage reduc- 














Fic. 2. Sketch illustrating prin- 
ciple of forced evaporation. 


ing valve, to which is connected a supply hose directly 
from the liquid oxygen supply at the-same pressure 
(min. 45 Ibs. to max. 80 Ibs. per sq.in. absolute) as it 
would be from a compressed oxygen cylinder after going 
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through the first-stage reducing valve on the B-L-B 
installations. 

Liquid oxygen (a) is in a galvanized sheet-iron con- 
tainer. A vacuum bottle can be used if the cost is 
justified, but because of the low temperature of the 
surrounding, perforated insulation container filled with 
glass wool (b) it is not a necessity. A similar insulating 
effect is produced by allowing cold gas to seep through 
the space (b) filled with glass wool. 

As the pressure rises simultaneously in the containers 
(a), (b), and (c), 45 Ibs. absolute is reached immedi- 
ately and the supply of gas is ready for use. If at any 
time, the pressure in (a), (b) and (c) goes above 80 Ibs. 
absolute, or any other predetermined pressure, the two 
safety valves (f) would release the excess pressure. 

In the meantime, the temperature of the container (a) 
has dropped and the rate of evaporation is lowered with 
consequent drop in the pressure in the containers (a), 
(b), and (c) as the result of continuous consumption. 

The moment this pressure approaches the required 
minimum of 45 Ibs., the valve (d) starts to close, raising 
slightly the pressure in the container (a). This pressure 
instantly forces the liquid oxygen in the tube (e) spilling 
liquid into the container (c) and producing a large 
volume of oxygen gas until enough pressure is produced 
to open valve (d). If at any time the pressure rises 
above the critical design pressure, safety valves (f) 
release the excess pressure. This process is continuous 
until all the liquid oxygen is used up. 
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Fic. 3. Schematic working assembly of liquid 
oxygen supply apparatus connected to B-L-B mask 
and second stage reducing valve. 


Fig. 3 shows the schematic working assembly of the 
apparatus as shown on Figs. 4 and 5. The additional 
items are: 


1. A perforated metal shield so that cold gas will 
uniformly surround the container (a). 
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2. A neck so arranged that liquid oxygen can be 
poured in from the outside without disturbing 
the setup. 

3. A visual gage to indicate the amount of oxygen in 
the container. The apparatus can be inclined 
or inverted. 

4. Safety valves. 

5. Automatic absolute pressure valve to restrict the 
gas flow from the container and spill the liquid 
in the tank. 

6. Gas filter. 

Starting by-pass valve. 


“J 


DESCRIPTION OF APPARATUS 


Gas Pressure Container 


The outside container is made of black iron, 0.065 
inch thick, but on actual installation it could be made 
of steel armor plate, one inch thick, thus utilizing the 
same material for armor and for the pressure container. 
The necessary heat conductivity to vaporize the liquid 
oxygen would not be impaired. 


Insulation 


The insulation material (glass wool) around the inside 
liquid reservoir is enclosed in a perforated, galvanized 
iron container and a low temperature is provided by 
the cold oxygen gas seeping through this insulation. 
Such insulation is not as good as a metal vacuum bottle, 
but, since a certain amount of steady flow of gas is 
needed, the poorer insulation acts favorably in this 
installation. Substitution of a metal vacuum bottle 
for the glass wool insulation is permissible and would be 
even better if a very small oxygen supply is needed. A 
vacuum bottle insulation inside of this apparatus would 
reduce the evaporation, and consequently the gas 
supply, to an approximate minimum value of 5 per cent 
per 24 hours as against approximately 10 per cent per 
hour for glass wool on small installations. A vacuum 
bottle would help if minimum oxygen volume is of 
major importance, but it also could be used to produce 
larger amounts of oxygen the same as with glass wool 
insulation. It is very questionable if, for all installa- 
tions, a vacuum bottle would be justified if we keep in 
mind that such vacuum bottles would have to be de- 
signed with no carbon absorption pads and that they 
must withstand pressure from both sides on the vacuum 
space. 


Liquid Oxygen Container 


The liquid oxygen container can be made out of gal- 
vanized steel but never of glass because of the risk of 
cracking during filling. This container is provided 
with a low-pressure relief valve in order not to allow a 
big difference of pressure in the inside of the container 
which would cause excessive spilling of oxygen and 
wastage of gas through the safety valve. The liquid 





Fic. 4. Experimental apparatus con- 


nected to B-L-B mask. 





Fic. 5. Experimental apparatus show- 
ing slight frost after evaporation of 6 liters 
of oxygen. 


container and glass wool insulation could be replaced 
by a vacuum bottle in one unit. 


Container Pressure Control Valve 


The container pressure control valve should be so 
constructed that it will open and close the gas outlet 
from the liquid container in order to spill liquid oxygen 
through the liquid outlet into the outside container very 
delicately, to maintain the pressure in the tank within 
the predetermined limits. Those limits might be set 
at will from one atmosphere to any desired pressure. 
This pressure is determined by the requirements of the 
second-stage reducing valve, passing gas to the mask. 
With further development of this apparatus there is a 
possibility of eliminating the second-stage reducing 
valve and substituting in its place a plain needle valve. 
This automatic inside pressure control valve is a very 
essential part of the apparatus and should provide a 
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smooth supply of gas and avoid excessive waste of oxy- 
gen through the safety valves. Two types of valves 
have been investigated, namely, a sylphon-operated 
sleeve valve (Fig. 6) and a spring-operated diaphragm 
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Fic. 6. Plunger valve. Absolute pressure control in 
oxygen container. 


valve. Both have advantages and disadvantages and 
detailed improvements in either one might make them 
of equal value. 


Containers for Liquid Oxygen Storage 


Containers for liquid oxygen storage are on the market 
in sizes from 5 liters to 100 liters with a rate of evapora- 
tion as low as 2 per cent in 24 hours (Figs. 7 and 8). 





Fic. 7. Small storage con- 
tainer. 


Handling of such liquid containers is not any more 
burdensome in supply bases than handling of gas cylin- 
ders since the weight of the liquid container is only a 
small per cent of the weight of the oxygen, while in gas 
cylinders the weight of the cylinders is a major part com- 
pared with the weight of the actual gas. It is of interest 
to note that the supply of oxygen gas in cylinders is 


always secured from a base where first liquid oxygen is 
delivered and then the gas cylinders are filled with 
oxygen gas by a process of evaporation and compression. 
In large scale operations using gas cylinders, liquid 
oxygen handling and supplying is unavoidable. 


Liquid Oxygen Producing Plants 


As stated above, the supply of cylinders of gaseous 
oxygen requires bases where liquid oxygen can be pro- 
duced or to which it can be shipped from central pro- 
ducing plants in special oxygen tank cars or in large size 
containers and then evaporated, compressed and filled 
in oxygen gas cylinders. 





; Storage Container 100 liter capacity ; 
Fic. 8. 


There is a belief that the supplying of liquid oxygen 
would be burdensome and almost an impossibility for 
such places as Army or Navy bases. Such bases would 
have to be supplied with a filling plant to which liquid 
oxygen would have to be delivered, and provision would 
have to be made to convert liquid into compressed gas 
in cylinders. To the surprise of the author, manufac- 
turers of oxygen equipment, on inquiry, submitted data 
showing that an oxygen plant, portable on a truck, could 
be easily manufactured. For example, a plant pro- 
ducing 5 liters of liquid oxygen per hour would be port- 
able on one truck and could be run with the same power 
plant as the truck itself or with an independent motor. 


Safety 


There have been a sufficient number of tests made 
which show that incendiary bullets shot through a liquid 
oxygen container, if the container is not manufactured 
with carbon absorption pads, will not present a fire 
hazard. As far as burning of surrounding objects is 
concerned, it makes no difference whether the oxygen 
gas comes from liquid or compressed oxygen gas. 
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Effect of Compass Location in Aircraft 
on Dynamic Errors 


WALTER WRIGLEY 
Sperry Gyroscope Company, Inc. 


SUMMARY 


A compass may be located, either for convenience or necessity, 
at some distance from the center of gravity of large modern air- 
craft. Such a location may give rise to dynamic errors in the 
readings of the compass due to the accelerations of roll, pitch 
and yaw of the airplane about its center of gravity. The accelera- 
tion of yaw is probably the most serious factor in producing such 
dynamic errors regardless of whether the compass is located in 
the nose, the tail or in a wing. 

The magnitude of the dynamic errors in compass readings in- 
creases in an approximately linear fashion as the distance of the 
compass from the center of gravity of the airplane is increased. 
For any particular location of the compass there will be a certain 
heading of the airplane for which the compass reading will have a 
maximum error. The dynamic errors may exist as magnitude 
effects of up to several degrees, or as phase effects of up to almost 
one hundred eighty degrees. 


INTRODUCTION 


o'er MAGNETIC compass, whose duty it is to furnish 
indications of azimuth to the pilot, may be located 
at some distance from the center of gravity of large 
modern aircraft. If, for reasons of convenience, it is 
located on the instrument panel in the pilot’s cockpit, 
the compass will generally be well forward of the air- 
plane’s center of gravity. On the other hand, it may 
be desirable to mount the compass as a remote reading 
instrument near the tail of the airplane or out in one 
of the wings in order to remove the effects of deviation 
due to the presence of magnetic materials or stray fields 
near the cockpit. Even if it were desired to mount the 
compass at the center of gravity of the airplane, such a 
location would be highly impractical since the center of 
gravity might be inaccessibly situated inside the struc- 
ture or in a space required for other purposes. In 
small craft the compass would probably never be lo- 
cated more than a few feet from the center of gravity, 
but in large ships it is more than likely to be mounted 
sufficiently far from the center of gravity to cause its 
location to have a marked effect on its performance. 

A magnetic compass indicates azimuth by showing 
the direction of the horizontal component of the earth’s 
magnetic field. Actually, the compass indicates on 
the compass card the direction of the projection of the 
earth’s magnetic field and for this reason the card must 
be pendulous. Except for dynamic effects the card is 
perpendicular to the direction of the resultant accelera- 
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tion acting on it and is, therefore, horizontal under the 
influence of gravity alone. Practical compasses now 
in use generally have such short natural periods of os- 
cillation about a horizontal axis in relation to the time 
required for the resultant acceleration of the instrument 
to change in direction that the compass card remains 
perpendicular to the resultant acceleration even for 
dynamic effects.' The card then tilts under the in- 
fluence of horizontal acceleration components. A tilt 
component of the compass card about the magnetic 
north-south axis, due to an east-west acceleration, 
introduces a projection of the vertical component of 
the earth’s magnetic field in the plane of the card. 


ue 
z CARD 


» NORMAL 






s HORIZONTAL 
PLANE 


Fic. 1. Projection of earth’s field in plane of 
card. Airplane banked 30° and headed north.! 


Fig. 1 shows how the projection H,, in the plane of the 
tilted compass card, due to the vertical component 
H, of the earth’s magnetic field H, combines with the 
horizontal component H), of the earth’s field to form a 
resultant field H, in the card. Starling? and Draper, 
Cook and McKay! have shown that the change in di- 
rection from the magnetic north of this resultant field 
in the compass card is the cause of dynamic errors in 
compass indications. 


ACCELERATIONS ACTING ON A Compass CARD 


When a compass is located at the center of gravity 
of an airplane, the instrument experiences only the ac- 
celerations of the airplane itself while in flight. Such 
accelerations accompany changes in linear velocity or 
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changes in the direction of flight during turns. When 
the compass is located away from the center of gravity 
of the airplane, the instrument will experience, in addi- 
tion to the above-mentioned accelerations, accelera- 
tions due to motions of the airplane about its center of 
gravity such as rolling, pitching or yawing. Fig. 2 
illustrates the origin of these accelerations experienced 
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IN TAIL 


AXIS OF ROLL ae 
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AXIS OF PITCH. 


Fic. 2. Various locations of a compass in an airplane. 
by a compass when mounted fore, aft or to one side of 
the center of gravity of the airplane. When mounted 
in position A forward of the center of gravity, or aft in 
position B, the compass is subject mainly to effects 
caused by yawing of the airplane about its vertical 
axis. The effect of pitching will add to or subtract 
from that of gravity and so will be of secondary impor- 
tance in tilting the compass. The effect of rolling will 
be small as long as the compass is not located more than 
a few feet from the axis of roll. When mounted in 
position C to one side of the airplane’s center of gravity, 
the compass is also subject mainly to yaw. The effect 
of rolling will add to or subtract from that of gravity, 
while the effect of pitching will be small as long as the 
compass is not located far from the axis of pitch. 

Horizontal acceleration due to yaw is thus seen to be 
the most important factor in tilting a compass that is 
located away from the center of gravity of an airplane, 
although rolling and pitching may also be considered 
contributing factors. Horizontal acceleration due to 
yaw is a tangential linear acceleration proportional to 
the product of the angular acceleration of yaw (a func- 
tion of the motion of an airplane) ‘and the distance the 
compass is located from the center of gravity. It is 
the latter factor that is of interest in this paper. 

The tilt of a compass card is dependent on the re- 
sultant of all the acceleration components—both the 
horizontal components which depend on the location 
of the compass in the airplane and the vertical com- 
ponents which are due to the force of gravity. 


COMPASS PERFORMANCE DURING A TURN 


It has been shown elsewhere by the present writer® 
that from a practical point of view the only horizontal 
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acceleration component effective in tilting a compass 
card during a turn is the centripetal acceleration of 
turning; the effects of accelerations due to yaw and 
roll are negligible. Since the centripetal acceleration 
acts at the center of gravity of an airplane, it follows 
that during a turn about the vertical the dynamic errors 
of a magnetic compass are substantially unaffected by 
the location of the compass in the airplane. 


COMPASS PERFORMANCE IN AN OSCILLATION ABOUT AN 
AXIS THROUGH THE CENTER OF GRAVITY OF AN AIR- 
PLANE 


To illustrate the effect of an oscillation on the dy- 
namic errors of a compass a sinusoidal yawing of the 
airplane about its vertical axis is chosen. This yawing 
motion produces a sinusoidally varying horizontal ac- 
celeration component which in turn causes a sinusoidal 
variation, about a horizontal axis, of the direction of 
the resultant acceleration acting on the compass. The 
compass card then executes a sinusoidal tilt about the 
same horizontal axis, and the north-south component 
of the tilt introduces a new magnetic field component 
in the plane of the compass card, tending to cause an 
error in the compass indications. The magnitude of 
the dynamic errors of the compass will be linearly pro- 
portional to the distance of the instrument from the 
center of gravity. This is one effect of compass loca- 
tion on dynamic errors. A second effect depends on 
the heading of the airplane. If the airplane in Fig. 2 
is flying in a general northerly or southerly direction, 
a compass located at points A or B, 1.e., fore or aft of 
the airplane’s center of gravity, will tilt about the 
north-south axis and so have dynamic errors due to 
location. However, if the compass is located in a wing 
at point C, the instrument will tilt about an east-west 
axis which does not lead to dynamic errors. But if the 
airplane is flying east or west, a compass located at 
point C will have dynamic errors, while one located at 
points A or B will not. For intermediate headings all 
compass locations will experience some north-south 
tilt components and hence lead to some dynamic 
errors. 

Fig. 3 shows the effect of location on the direction of 
the resultant magnetic field in the plane of the compass 
card during a sinusoidal oscillation of the airplane. 
In Fig. 3 the direction of the resultant field is plotted 
against the time for one period of yaw. It is assumed 
that the period of yaw is four seconds, the amplitude of 
yaw is ten degrees, and the airplane is flying north 
over New York. It is seen in the figure that the re- 
sultant field points north when the compass is located 
at the center of gravity of the airplane, but the field 
oscillates about north in a period equal to the period of 
yaw for non-center-of-gravity locations. For the 
particular case cited the maximum variation in field di- 
rection is about ten degrees when the compass is located 
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TIME IN SECONDS 


COMPASS MOUNTED 20FT. FORWARD OF CG 
———— COMPASS MOUNTED | OFT.FORWARD OF CG 
—-—--COMPASS MOUNTED AT CG 

PERIOD OF AIRPLANE YAW 4 SECONDS 
AMPLITUDE OF AIRPLANE YAW 10° 
AIRPLANE FLYING NORTH OVER NEW YORK 
NOTE: FIELD HEADING WILL BE 180° OUT 
OF PHASE WHEN THE COMPASS IS 
LOCATED AFT OF THE CG 





Fic. 3. Direction of resultant magnetic 
field in the plane of the compass card as a 
function of the location of the compass in the 
airplane during a sinusoidal oscillation in yaw. 


ten feet forward of the center of gravity and about 
nineteen degrees when located twenty feet forward. 

The response of a compass to this field oscillation is 
shown in Fig. 4. The reading of the compass is plotted 
against the time for one period of yaw. The compass 
used is an average of several standard instruments and 
has an undamped natural period about its axis of sym- 
metry of ten seconds;! and, in addition, is critically 
damped (aperiodic). As seen in Fig. 4, the compass 
gives correct indications of the heading of the airplane 
when the instrument is located at the center of gravity 
of the airplane. When located away from the center 
of gravity, the compass readings are in error in both 
magnitude and phase. For the particular case chosen 
the compass readings show about seventy-five per cent 
of the true magnitude, but are one hundred thirty-five 
degrees out in phase when the compass is located for- 
ward of the center of gravity. The readings are about 
double in magnitude, but only forty-five degrees out in 
phase when the compass is located abaft the center of 
gravity. When mounted in the wing, the compass 
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Fic. 4. Dynamic errors in compass indications 
as a function of location of the compass in an air- 
plane during a sinusoidal oscillation in yaw. 


gives correct indications. It should be noted here that 
the curves of Fig. 4 are based on the assumption of 
steady state oscillations. 

If the airplane were flying south under the conditions 
given in Figs. 3 and 4, the effect of compass location 
would be reversed. In other words, the compass indi- 
cations are about double in magnitude and almost cor- 
rect in phase when the instrument is mounted in the 
nose. A compass located in the right wing will perform 
for a westerly course as it does for the northerly course 
as illustrated in Fig. 4 when mounted in the nose. 
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WILBuR Wricut MemoriAL LECTURE 


Juan T. Trippe, M.I.Ae.S., President of Pan American Airways 
System, and a member of the Advisory Board of the Institute, 
delivered the Twenty-Ninth Wilbur Wright Memorial Lecture 
before the Royal Aeronautical Society in London on June 17, 
1941. 

Mr. Trippe’s lecture on ‘Ocean Air Transport’’ described the 
history and development of Pan American’s transoceanic service, 
the technical methods used in planning and controlling over-ocean 
flights, and the research on equipment and methods necessary for 
their improvement. 

The invitation to deliver the annual Wilbur Wright Memorial 
Lecture is one of highest honors conferred by the Royal Aero- 
nautical Society. The lecturers chosen each year are alternately 
British and foreign aeronautical specialists. 


AWARD FOR Foc Stupy 


The paper ‘“‘Fog Study of Miami, Florida,’’ by J. R. DeHart, 
published in this issue of the Journal, is a report on the research 
which won for its author an award from the Air Transport 
Association of America. 

Coincident with the reading of his paper at the meeting of the 
American Meteorological Society in Miami on May 8, it was an- 
nounced that Mr. DeHart had been selected, for this study, as 
the first recipient of this annual award with a cash prize of $500. 
It is conferred by the Air Transport Association to employees of 
airlines in the United States ‘‘for the best original research having 
practical application to meteorology and dispatching.’”’ Recipi- 
ents are selected by the Operations Committee of the Association, 
the members of which are the vice-presidents in charge of opera- 
tions of the nation’s airlines. 

Mr. DeHart, a Technical Member of the Institute, is a meteor- 
ologist on the staff of Eastern Air Lines stationed at Miami, 
Florida. His paper was part of a general research program 
started by Eastern Air Lines last summer under the supervision 
of Joseph J. George, Chief of its Meteorological Department. 

DeHart, who is 26, was graduated from the University of Mi- 
ami with a Bachelor of Science degree in 1937. He then took a 
year of post-graduate work in electrical engineering at Georgia 
School of Technology followed by a course in meterology at the 
Boeing School of Aeronautics. He has been with Eastern Air 
Lines since the spring of 1939. 


New CorporaTE MEMBERS 


The “‘Appreciation Pages’’ in the front of the Journal now list 
six additional Corporate Member companies since the last an- 
nouncement of new Corporate Membets made in the March 
issue. 

These companies affiliated with the Institute since that date 
are: 


Associated Aviation Underwriters, New York, N. Y., aeronau- 
tical insurance. 

Eaton Manufacturing Company, Cleveland, Ohio, manufac- 
turers of aircraft engine parts and accessories, pumps, valves 
and actuating cylinders. 

Casey Jones School of Aeronautics, Newark, N. J., and its 
affiliated Academy of Aeronautics, La Guardia Field, N. Y., 
which give courses in aeronautical engineering, master me- 
chanics’ training, airplane design and construction, airplane 
mechanics and maintenance. 
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North American Aviation, Inc., Inglewood, Calif., and Dallas, 
Texas, constructors of the military aircraft bearing its name. 

Packard Motor Car Company, Detroit, Michigan, which is start- 
ing the production of Rolls-Royce Merlin XX aircraft engines. 

Vultee Aircraft, Inc., Vultee Field, Calif., and Nashville, Tenn., 
with its Stinson Aircraft Division, Wayne, Mich., manufac- 
turing Vultee military airplanes and Stinson commercial air- 
planes. 


A list of the products and services of Corporate Member com- 
panies is now published in expanded form in the Directory in the 
‘Aeronautical Review Section” of the Journal. It will be ap- 
preciated if Institute members and readers of the Journal and 
“Aeronautical Review’’ will give consideration to the items 
listed in this Directory whenever possible. 

Corporate Membership was established in 1936 for companies 
which wished to participate in the activities of the Institute and 
contribute to the publication of the Journal. It is limited to 
“manufacturers or operators of aircraft or producers of equip- 
ment, services or supplies used in aeronautics.’’ The membership 
dues contributed by Corporate Members have enabled the Insti- 
tute to expand its services over the past years and help make 
possible the publication of its technical Journal. 

The greater contribution to the progress of aeronautics made 
possible by their support of the Institute is appreciated by all 
members. 


Honorary DEGREES TO INSTITUTE OFFICERS 


Major James H. Doolittle, F.I.Ae.S., immediate Past President 
of the Institute, and Ralph S. Damon, A.F.I.Ae.S., member of the 
Institute’s Council, were the recipients of honorary degrees last 
month. 

Mr. Damon, President of Republic Aviation Corporation and 
recently Vice-President in charge of Operations for American Air- 
lines, received the honorary degree of Doctor of Engineering 
from Clarkson College of Technology at its forty-second com- 
mencement exercises on June 2, in Potsdam, N. Y. The degree 
was conferred “for meritorious contributions and outstanding 
advances in the fields of engineering and science.”’ 

The honorary degree of Doctor of Engineering was conferred 
upon Major Doolittle, who is serving in the Procurement Division 
of the U.S. Army Air Corps as coordinator of aircraft production 
in the automotive industry, by the Polytechnic Institute of 
Brooklyn. At the school’s annual Corporation dinner, preceding 
commencement exercises, Major Doolittle was presented for the 
degree by Lester D. Gardner, Executive Vice-President of the 
Institute of the Aeronautical Sciences. In his presentation 
speech, Major Gardner read letters of congratulations from Presi- 
dent Robert G. Sproul of the University of California and 
President Karl T. Compton of the Massachusetts Institute of 
Technology, at which schools Major Doolittle had received the 
degrees of A.B. and D.Sc., respectively; from Harry F. Guggen- 
heim, who recalled the pioneering work in blind landing research 
done by Major Doolittle for the Guggenheim Foundation; from 
Alexander Fraser, President of the Shell Oil Company for which 
Major Doolittle had until recently been Aviation Director; and 
from Major General H. H. Arnold, as Chief of the U.S. Army Air 
Corps, in which he referred to Major Doolittle’s record in the Air 
Corps before his resignation in 1931 and expressed pleasure at 
again having his services available in its present procurement 
program. 

At the commencement exercises on June 11, President H. S. 
Rogers of Polytechnic Institute of Brooklyn read the following 
citation: 
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‘‘Major James Harold Doolittle, speed and safety of flight have 
focused your energies and challenged your genius. You, above 
all others, have contributed to their development in commerce 
and in defense. You have conquered distance and defied time. 
Your accomplishments have been made, not through good fortune 
befalling dauntless courage, but through careful use of advanced 
scientific knowledge and of precise engineering measurements. 
Your flights and experiments leading up to them have made 
aviation history and have kept this nation in the forefront of 
aviation. You have brought to this exacting field of endeavor 
the highest achievements of intellect and the finest qualities of 
courageous spirit. 

“Acclaiming your pre-eminence in aviation and applauding 
the manner in which you carry the signal honors you have re- 
ceived, in the name of Polytechnic Institute I confer upon you the 
degree, Doctor of Engineering, with all its rights and privileges.” 


News oF INsTITUTE MEMBERS 


Cecil W. Armstrong, M.I.Ae.S., Assistant Professor Engineer- 
ing, University of California, is employed for the summer months 
as a Senior Research Engineer at the Lockheed Aircraft Corp. 

Ensign Dan L. Belvin, U.S.N.R., Technical Member, who was 
formerly assigned as an Engineering Officer in the Maintenance 
Division, U.S. Naval Air Station, Pensacola, is now on duty in 
Washington at the Bureau of Aeronautics, Navy Department. 

Major General George H. Brett, who, at the time he was elected 
an Honorary Member of the Institute, had been serving as Acting 
Chief of the U.S. Army Air Corps, was appointed on May 23 as 
Chief of the Air Corps for a regular four year term. 

Ray W. Brown, M.I.Ae.S., has left his position as General 
Manager of the Aviation Department of The General Tire & 
Rubber Company to go on active duty in the Bureau of Aero- 
nautics, Navy Department. 

William A. M. Burden, M.I.Ae.S. and a Vice-President of the 
Institute, has resigned as a Director of National Aviation Corp. 
He is connected with the recently established unit of the Recon- 
struction Finance Corporation concerned with financial aid to 
Latin American aviation. 

John E. Cameron, Technical Member, has left the employ of 
Boston-Maine Airways to join the engineering staff of Boeing Air- 
craft Company as a layout man. 

Lt. Clinton D. Case, U.S.N., A.F.I.Ae.S., has been detailed to 
duties aboard the U.S.S. Albemarle. He had previously been 
Officer-in-Charge, Power Plant Division, Ground School, U.S. 
Naval Air Station at Pensacola. 

Charles Healy Day, A.F.I.Ae.S., has recently been appointed 
to the Canadian Government’s Department of Munitions and 
Supply as Supervisor of Production, Aircraft Branch, with head- 
quarters in Ottawa. Mr. Day established an aircraft factory for 
the Chinese Government at Shiuchow in 1934 and for the next 
six years was technical adviser on aircraft production to the 
National Government of China. 

Dr. J. P. Den Hartog, F.I.Ae.S., Associate Professor of Ap- 
plied Mechanics, Harvard University, has been assigned to active 
duty in Washington as a Lieutenant Commander, U.S. Naval 
Reserve, at the Bureau of Ships, Navy Department. 

George A. Doole, Jr., M.I.Ae.S., has been appointed Opera- 
tions Manager for Aerovias Nacionales de Colombia of the Pan 
American Airways System at Barranquilla. He was formerly 
Chief Pilot of AVIANCA. 

At the annual meeting of the Aviation Writers Association 
held in Louisville, Ky., June 1-4, 1941, Michael Froelich, M.I.- 
Ae.S., News Editor of Aero Digest was elected Treasurer and 
Carl Norcross, M.I.Ae.S., Managing Editor of Aviation, Secre- 
tary of the Association. 

L. W. Glowa, Technical Member, formerly with the Sales 
Department of American Airlines in New York, is now an Air- 
ways Engineer for the U.S. Army Air Corps in Washington. 
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A. I. Lodwick, M.I.Ae.S., formerly with Aviation Manufac- 
turing Corporation, has recently established the Lakeland School 
of Aeronautics at Lakeland, Florida, in partnership with H. S. 
Darr. 

R. C. McKenzie, Technical Member, formerly with De Havil- 
land Aircraft of Canada Ltd. is now serving as an aeronautical 
engineer with the Royal Canadian Air Force at Montreal. 

Oliver L. Parks, M.I.Ae.S., President of Parks Air College, was 
recently elected President of the Southwest Aviation Conference. 
L. G. Fritz, M.I.Ae.S., Vice-President in charge of Operations for 
TWA, was elected Vice-President of the Conference. 

E. Robert Reff, M.I.Ae.S., formerly Structural Development 
Engineer with the Airplane Division of the Curtiss-Wright 
Corporation at Buffalo, has joined the engineering staff of 
Northrop Aircraft, Inc. 

Laurance S. Rockefeller, M.I.Ae.S., has been elected a member 
of the Executive Committee of the Board of Directors of Mc- 
Donnell Aircraft Corp. 

A. P. Taliaferro, Jr., M.I.Ae.S., has resigned from the vice- 
presidency of Aeronautical Securities, Inc., to go on active service 
in the U.S. Navy as Lieutenant Commander. 

At a recent meeting of the Board of Directors of Wright 
Aeronautical Corp., Philip B. Taylor, F.I.Ae.S., Assistant General 
Manager, was elected Vice-President, and William D. Kennedy, 
M.I.Ae.S., Sales and Service Manager, was elected Vice-President 
in charge of Sales. 

Herbert V. Thaden, A.F.I.Ae.S., formerly Assistant Sales 
Manager of Hayes Manufacturing Corp., has been appointed 
Manager of the Commercial Department of Duramold Aircraft 
Corp., New York. 

Charles F. Thomas, M.I.Ae.S., Technical Sales Engineer with 
Lockheed Aircraft Corp., is now in charge of the company’s 
office in Dayton. 

Lt. Comdr. D. W. Tomlinson, M.I.Ae.S., has been appointed 
Commanding Officer of the U.S. Naval Reserve Aviation Base at 
Fairfax Airport, Kansas City. He is on indefinite leave of ab- 
sence from his office as Vice-President in charge of Engineering 
for TWA, Inc. 

Howard M. Willson, Technical Member, has been transferred 
from the Cleveland Airport office of American Airlines to the 
company’s office at Lockheed Air Terminal, Burbank, California. 


Honorary MEMBERS ELECTED 


Robert A. Lovett, Assistant Secretary of War for Air, and 
Merrill C. Meigs, Chief of the Aircraft Section, Office of Produc- 
tion Management have been elected to Honorary Membership 
in the Institute. 

Merrill Church Meigs was born in Malcolm, Iowa, November 
25, 1883. He studied at Iowa Business College and in business 
courses at the University of Chicago. For the past thirty years 
he has been prominent in advertising and publishing activities in 
Chicago. He was successively publisher of the Chicago Herald 
and Examiner from 1926 to 1929; of the American Weekly from 
1930 to 1933 and of the Chicago American from 1933. Since his 
first airplane flight in 1928 he has been an enthusiastic amateur 
pilot, having obtained a Transport Pilot’s license and accumulated 
close to 1300 hours of solo flying time. He has been chairman of 
the Chicago Aeronautics Commission and Illinois Aeronautics 
Commission. He was named to head the Aircraft Section of 
OPM in January of this year. 

Robert Abercrombie Lovett was born September 14, 1896, 
at Huntsville, Texas, and received a B.S. degree from Yale 
University in 1918 followed by studies in law and business ad- 
ministration at Harvard University. During the World War he 
was a pilot in the U.S. Naval Air Service, rose to the rank of 
Lieutenant Commander and was awarded the Navy Cross for 
his service overseas. His extensive experience in finance and 
industry began in 1921 with a clerkship in the National Bank of 
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Commerce in New York and has included membership on the 
boards of directors of several railroads and manufacturing com- 
panies, insurance and banking institutions. He has also been 
active in numerous charitable organizations in New York City. 
Mr. Lovett was a partner in the firm of Brown Brothers, Hari- 
man & Co. until called in December, 1940, to serve as a special 
assistant to Secretary of War Stimson. His accomplishments in 
this position in handling executive matters pertaining to the ex- 
pansion of the Air Corps led to his appointment by President 
Roosevelt on April 10, 1941, to the long vacant post of Assistant 
Secretary of War for Air. In this office his efforts will be of great 
importance in the unification of authority within the Air Corps 
and better coordination with the entire national defense effort. 


].Az.S. StuDENT BRANCH AWARDS 


Recipients of the two Student Branch awards, the I.Ae.S. 
Student Branch Scholastic Award and the I.Ae.S. Student Branch 
Lecture Award, in addition to those listed in the June issue of 
the Journal, are as follows: 

Student Branch Scholastic Award: B. Shreve Hutchins, 
Agricultural and Mechanical College of Texas; Frank T. Lewand, 
University of Detroit; Robert G. Urquhart, Harvard University; 
Leon R. Wosika, University of Michigan; Robert W. Bratt, 
University of Minnesota; Walter Walling, Oregon State College; 
LeRoy Krouse, Tri-State College; and Victor M. Ganzer, Uni- 
versity of Washington. 

Student Branch Lecture Award: Jack Holt, Agricultural and 
Mechanical College of Texas; William S. Holmes, Polytechnic 
Institute of Brooklyn; Frank T. Lewand and Frank J. Bajkowski 
(jointly), University of Detroit; J. Norman Ball, Jr., Harvard 
University; Leslie J. Trigg, University of Michigan; James A. 
Butcher, University of Minnesota; Raymond Pederson, Oregon 
State College; Warren LeRoy, Tri-State College; and William 
F. Kaseburg, University of Washington. 


PERSONNEL OPPORTUNITIES 


The Personnel Bureau serves individual members, as well as 
organizations seeking to employ aeronautical specialists. Any 
member or organization may have requirements listed without 
charge. 


Wanted 


Foreman—Balloons—completely familiar with construction 
of lighter-than-air craft. Location in East. State experience and 
salary desired. Address reply to Box 135, Institute of the Aero- 
nautical Sciences. 


STUDENT BRANCHES 


Boeing School of Aeronautics. At a recent meeting, the 
following new officers were elected: Andrew Smisko, Chairman; 
Charles Strausser, Vice-Chairman; and Claude Williams, Secre- 
tary-Treasurer. , 

University of California. Charles W. Harper, Chairman, ap- 
pointed Bradford Wick and Arthur Hughes as a committee of two 
to organize the Student Branch for the Fall semester. 

University of Detroit. George H. Tweney, Instructor in 
Aeronautics, presented the I.Ae.S. Branch Lecture on ‘‘Stainless 
Steel and its Aircraft Uses” by Dr. Michael Watter, at a recent 
meeting. New officers were elected as follows: James Hailstones, 
Chairman; James Rice, Vice-Chairman; Wilfred Bush, Treas- 
urer; and Louis B. Peltier, Secretary. Prof. George J. Higgins 
was re-elected Faculty Sponsor. 
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University of Michigan. The following men have been elected 
to serve as officers for the school year 1941-42: Henry T. Field- 
ing, Chairman; J. Schwarzbach, Treasurer; and Arthur Billet, 
Secretary. Elections for the office of Vice-Chairman and Honor- 
ary Chairman will take place in the Fall. 

New York University. At a meeting held on May 7, officers 
were elected as follows: James Smith, Chairman; Ralph Nun- 
ziato, Vice-Chairman; and Frank Shieflet, Secretary-Treasurer. 

Oregon State College. On May 12, officers were elected as 
follows: Lawrence Robertson, Chairman; Forbes McRae, Vice- 
Chairman; and Don Jacoby, Secretary-Treasurer. Prof. Ben- 
jamin F. Ruffner was re-elected Honorary Chairman. 

Ryan School of Aeronautics. A dinner meeting was held on 
April 18 at which time Walter B. Sayner read a paper on ‘“‘The 
Assisted Take-Off.”” Ata later meeting, held on May 23, W. T. 
Immenschuch gave a paper on “The Work of the Designing 
Engineer.” 

Virginia Polytechnic Institute. Officers for the next school 
year were elected as follows: J. W. Martin, Jr., Chairman; E. F. 
Burner, Vice-Chairman; A. Hertzberg, Secretary; and W. W. 
Lumpkins, Treasurer. 


CHANGES IN MEMBERSHIP 


The following changes in the membership of the Institute 
have occurred since the publication of the previous issue of the 
Journal. 


ELECTED TO MEMBER GRADE 


Dufort, Robert Harold, B.S. in Ae.E.; Head of Stress Dept., 
Goodyear Aircraft Corp. 

Markey, Richard Joseph, Final Assembly Production Planning, 
Vought-Sikorsky Aircraft, United Aircraft Corp. 

Montgomery, Wiley Pope, B.S. in Ae.E.; In charge of Structural 
Research, Vought-Sikorsky Aircraft, United Aircraft Corp. 

Shields, Bert A., Aviation Writer; Ground Flying Specialist and 
Lecturer, Columbia University; New York University. 

Stout, Ernest Gordon, Ae.E.; Engineer in charge of Hydro- 
dynamics, Consolidated Aircraft Corp. 


ELECTED TO INDUSTRIAL MEMBER GRADE 


Favorite, Lewis P., Sales (Engineering Service), Aluminum Co. 
of America. 

Warner, Percy deForest, II, Asst. Chief, Ground School Sec., 
Standards Div., Civilian Pilot Training, C.A.A. 


ELECTED TO TECHNICAL MEMBER GRADE 


Allen, Harold Frederick, M.S. in E.E.; Instructor, University of 
Michigan. 

Brown, Chester Harvey, Jr., B.S.; Junior Aero. Engineer, Grum- 
man Aircraft Engineering Corp. 

Carder, Alden Bivins, B.S. in Ae.E.; Flight Test Engineer, 
Douglas Aircraft Co. 

Delany, Noel Keesling, B.S. in Ae.E.; Jr. Aero. Engineer, 
Ames Aero. Lab., N.A.C.A. 

Dunn, Louis Gerhardus, Ph. D.; Instructor in Aeronautics, Calif. 
Inst. of Technology. 

Fruauff, George Louis, Jr., Bench Metal Man, Brewster Aero- 
nautical Corp. 

Graham, William Henry, M.A.; Assoc. Editor, Aeronautical 
Review, Institute of the Aeronautical Sciences. 

Hook, Franklin Louis, B.S. in Ae.E.; Draftsman, The Glenn L. 
Martin Co. 

Johnson, Paul Alford, B.S.E.; Vibration & Stress Engineer, 
Boeing Aircraft Co. 

Katzenberger, Edward Frederick, B.M.E.; Senior Layout Engi- 
neer, Brewster Aeronautical Corp. 
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Kruchek, Keith Cornelius, B.S. in Ae.E.; Engineer (Structures), 
Boeing Aircraft Co. 

Lindblom, Charles William, S.B. in Ae.E.; Shop Liaison Engi> 
neer, Lockheed Aircraft Corp. 

Markowitz, Jesse, M.S. in M.E.; Instructor, School of Tech- 
nology, City College of New York. 

Myers, Oscar Bruce, Jr.; B.S. in M.E.; Military Pilot, U.S. 
Army Air Corps. 

Nicol, William Franklin, Chief Engineer, Braniff Airways, Inc. 

Nietsch, Herman Erich, S.B.; Aero. Engineer, The Glenn L. 
Martin Co. 

Paradiso, Michael Antonio, Layout Design Engineer, Vought- 
Sikorsky Aircraft, United Aircraft Corp. 

Roberts, Phillip H., B.S.E.; Senior Detail Draftsman, Lockheed 
Aircraft Corp. 

Rohl, Henry Thomas, B.S. in Ae.E.; Jr. Aero. Engineer, U.S. 
Army Air Corps, Materiel Div., Wright Field. 

Sawyer, James Garnett, B.S. in M.E.; Airplane Division (St. 
Louis Plant), Curtiss-Wright Corp. 

Smith, George Campbell, B.S. in Ae.E.; Engineer Officer, Royal 
Canadian Air Force. Canada. 

Street, Robert Elliott, Ph.D.; Instructor in Math., Rensselaer 
Polytechnic Inst. 

Werner, Louis Fred, President & Gen. Manager, Models, Inc. 


Williams, Richard Samuel, Ae.E.; Technical Data Engineer, 
Intercontinent Aircraft Corp. 

Wright, Frederick Hamilton, Aerodynamics Dept., Douglas Air- 
craft Co. 


NECROLOGY 
HERBERT MARSH 


Herbert Marsh, a member of the Institute, Superintendent of 
Instruction at the Boeing School of Aeronautics, died on June 10, 
1941, at his home in San Leandro, California. 

Mr. Marsh was born in Lancashire, England, December 22, 
1886, and educated in the Viduct School of Technology there. 
At the outbreak of the World War in 1914 he left the Daimler 
Motor Company of Coventry, England, to become Chief Examiner 
in the Aeronautical Inspection Department of the British Air 
Ministry. With acommission as Captain, he was Chief Technical 
Officer in the Materiel Division (Research Supply) of the Royal 
Flying Corps and the Royal Air Force from 1916 to 1919. Follow- 
ing this he was a Works Inspector for the Air Ministry. 

On coming to the United States, he joined the National Parks 
Airways in 1929 as Assistant Manager of Operations. He became 
an Instructor on the staff of the Boeing School of Aeronautics in 
1935 and had been with the school from then until the time of his 
death. 





Book Review 


Theodor von Karman Anniversary Volume; Contributions to 
Applied Mechanics and Related Subjects, by the Friends of 
Theodore von K4arman on his Sixtieth Birthday; California 
Institute of Technology, Pasadena, 1941; 337 pages, $3.75. 

The twenty-six scientific papers presented in this volume are 
preceded by a short biographical note by Clark B. Millikan witha 
list of ninety-one important publications by von Karm4n in the 
fields of Applied Mathematics, Physics, Strength of Materials, 
Stress Analysis, Theory of Elasticity, Elastic Stability, Mono- 
coque Structures, Vibrations, Mechanics of Ideal, Viscous, and 
Compressible Fluids, Turbulence, Aerodynamics of Aircraft 
Hydrodynamics of Planing Surfaces, and Heat Transfer. 

The Table of Contents of this volume best indicates its char- 


acter: 


H. Bateman, California Institute of Technology: 
Some Definite Integrals Occurring in Aerodynamics 
Max M. Munk, Catholic University of America: 
On the Geometry of Streamlining 
J. C. Hunsaker, Massachusetts Institute of Technology: 
Dimensional Analysis and Similitude in Mechanics 
Theodore Theodorsen, National Advisory Committee for Aero- 
nautics: 
Impulse and Momentum in an Infinite Fluid 
Boris A. Bakhmeteff, Columbia University: 
Coriolis and the Energy Principle in Hydraulics 
H. U. Sverdrup, Scripps Institution of Oceanography of the Uni- 
versity of California: 
The Influence of Bottom Topography on Ocean Currents 
W. F. Durand, Stanford University: 
Ionization as a Factor in Fluid Mechanics 
Hugh L. Dryden, National Bureau of Standards: 
Isotropic Turbulence in Theory and Experiment 
J. L. Synge and W. Z. Chien, University of Toronto: 
The Intrinsic Theory of Elastic Shells and Plates 
Francis D. Murnaghan, Johns Hopkins University: 
The Compressibility of Solids under Extreme Pressures 
F. Zwicky, California Institute of Technology: 
Hydrodynamics and the Structure of Stellar Systems 


H. M. Westergaard, Harvard University: 
On the Elastic Distortion of a Cylindrical Hole by a Localized 
Hydrostatic Pressure 
Paul S. Epstein, California Institute of Technology: 
On the Absorption of Sound Waves in Suspensions and Emul- 
sions 
R. Courant, New York University: 
On a Method for the Solution of Boundary- Value Problems 
Karl Arnstein, Goodyear Aircraft Corporation: 
The Engineering Treatment of Ring or Wheel Problems 
A. Einstein, V. Bargmann, and P. G. Bergmann, The Institute 
for Advanced Study: 
On the Five-Dimensional Representation of Gravitation and Elec- 
tricity 
S. Timoshenko, Stanford University: 
The Forced Vibrations of Tie-Rods 
Th. Troller, Daniel Guggenheim Airship Institute: 
A Note on the Effect of the Wind-Tunnel Size on Pitching 
Moments 
A. Nadai, Westinghouse Research Laboratories: 
The Creep of Metals under Various Stress Conditions 
K. O. Friedrichs, New York University: 
On the Minimum Buckling Load for Spherical Shells 
William Bollay, Harvard University: 
The Theory of Flow through Centrifugal Pumps 
Frank L. Wattendorf, Wright Field: 
The Ideal Performance of Curved-Lattice Fans 
L. H. Donnell, Illinois Institute of Technology: 
Stress Concentrations Due to Elliptical Discontinuities in Plates 
under Edge Forces 
Hans Reissner, Illinois Institute of Technology: 
On Lubrication Flow with Periodic Distribution between Pre- 
scribed Boundaries 
R. v. Mises, Harvard University: 
Some Remarks on the Laws of Turbulent Motion in Tubes 
Wolfgang B. Klemperer, Douglas Aircraft Company: 
Stress Pattern Crazing 
J. C. HUNSAKER 





SUGGESTIONS FOR CONTRIBUTORS TO THE JOURNAL 
of thee AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites 
both members and non-members from any country to 
submit papers for publication in the Journal of the 
Aeronautical Sciences. The Journal, following the 
practice of other scientific publications, does not pay for 
contributions. 


The following directions for the preparation of papers, 
if followed by authors, will save correspondence, avoid 
the return of papers for changes, minimize the work of 
preparation for the printer, and save the expense due 
to the charges made for ‘‘author’s corrections.” 





ManuscrirTs: Papers must be written in English, in original 
typewriting on one side only of white paper sheets, consecutively 
numbered, and be double or triple spaced with wide margins. 
Manuscripts should be prepared with great care so that they 
will be typographically accurate. Paragraphing should be given 
special attention. Papers should be written in the third person, 
reference to the writer being made as ‘“‘the author.” Avoid the 
use of the words “I,” “we,” and “you.” Blueprint copies of 
papers are unacceptable as it is impossible to mark directions to 
the printer on them. Correcting, changing, or adding to papers 
after they are in type is costly. It is, therefore, imperative that 
papers submitted be in final form. Typographical errors may 
be corrected on proofs, but if authors wish to add material, they 
may do so at their own expense. In mailing, drawings may be 
rolled, but manuscripts should be sent flat. Send by first class 
mail (register if you wish for your own protection) to the Secre- 
tary, Institute of the Aeronautical Sciences, 1505 RCA Bldg. 
West, Rockefeller Center, New York City. All manuscripts will 
be examined by the Editorial Committee and by the Editor. 
Authors will be advised as promptly as possible (usually two to 
three weeks) whether the paper is acceptable for publication. 


TitLEs: The title of the paper should be brief. The name and 
initials of the author should be written as he prefers. The use 
of the full name of an author is advocated because of the fre- 
quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
particularly important for large annual indexing and abstract- 
ing services. All titles and degrees or honors are omitted. The 
name of the organization with which the author is associated 
should be placed after his name on the same line. The date on 
which the paper is received will be inserted by the Editor. 


SUMMARIES OR ABSTRACTS: An abstract to be printed at the 
beginning should accompany each article. It should be ade- 
quate as an index and asa summary. It should contain a state- 
ment of major conclusions reached, since summaries in many 
cases constitute the only source of information used in compiling 
scientific reference indexes. Abstracts printed in other jour- 
nals, especially foreign, in most cases, consist of summaries 
from printed papers. Thesummary should explain as adequately 
as possible the major conclusions to a non-specialist in the 
subject. The summary should contain from 100 to 300 words, 
depending on the length of the paper. 


SuB-HEADINGS: Sub-headings should be inserted by the au- 
thor at frequent intervals. The work of editorial preparation 
will be simplified by the author providing many sub-headings. 
Owing to the breaking of columns and the insertion of illustra- 
tions, some of the sub-headings may have to be omitted. 


SHORTENING OF PAPERS: Some papers, at the end, fill in anly a 
portion of a page. This leaves much wasted blank space as 
succeeding articles are started at the top of a page. Authors 
should indicate by notation on the left-hand side of the page 
what matter may be omitted when “run overs” occur. This 
request is important as the Journal cannot afford in the future, 
as it has in its earlier issues, to have blank half pages or more 
at the end of papers. 


MATTER USUALLY DELETED: Acknowledgments of assistance 
in preparation of paper, except by collaborators. Photographs 
or illustrations of little technical interest and not showing ad- 
vances in general practice. Too detailed tabular matter (gen- 
eral results of such tables may be included in the text). Lengthy 
descriptions of materials or processes or of preliminary experi- 
ments or theories which preceded final results; salient features 
only are of interest. 


REFERENCES AND FOOTNOTES: References should appear as 
footnotes only, numbered consecutively, grouped together at the 
end of the manuscript. The arrangement should be as follows: 
(for books)—! Durand, W. F., Aerodynamic Theory, Vol. 1, p. 23; 


Julius Springer, Berlin, 1934. (For magazines)—Englund, 
C. R., Crawford, A. B., and Mumford, W. W., Some Results of a 
Study of Ultra-Short-Wave Transmission Phenomenon, Proc. I. 
R. E., Vol. 20, No. 12, pp. 481 and 482, March, 1933. Please 
give Author, Title, Volume, page, publisher and date of publi- 
cation as indicated. Omission of one required fact causes much 
extra editorial work and possible inaccuracies. All references 
are grouped at the end of the article. 

ILLUSTRATIONS: Illustrations should accompany manuscripts 
and each should always be referred to in the text, preferably by 
number. Drawings or graphs should not be larger than 12 X 16 
inches, and must be made with jet black India Ink on white 
paper or tracing cloth, the latter being preferred. Do not use 
typewriter for lettering. The smallest lettering on 8 X 10 inch 
figures should be no less than '/, inch high. Cross-section 
paper (white with black lines) may be used, but should not have 
more than 4 lines per inch. If finer ruled paper is used, the 
major division lines should be drawn in with black ink, omit- 
ting the finer divisions. In the case of finely ruled paper, only 
blue-lined paper can be accepted. Tracing paper and blue- 
prints are not acceptable. Lettering and all markings must 
be large enough to be readable after reduction. Mail rolled or 
flat, never fold. Drawings which cannot be reproduced (includ- 
ing pencil drawings) will be returned to the author for redraw- 
ing, thus delaying publication of the paper. Photographs 
should be very distinct and show clear black and white con- 
trasts. They must be on glossy white paper. Avoid round and 
oval photographs. 

CAPTIONS AND LEGENDS: Legends or captions must accom- 
pany each drawing or photograph submitted. If written on the 
drawing or photograph, they should be placed below and well out- 
side the part to be reproduced. It is better to place them on 
separate sheets of paper pasted to the back of the drawings or 
photographs. Each table should have a caption such as Table 
1, Table 2, Table 3, etc. Captions should be complete in them- 
selves so as to make the data intelligible to the reader without 
reference to the text. A duplicate list of captions for figures 
should be included as the last page of the manuscript. Use 
“Fig. 1” (not Figure 1), Figs. 3 and 4, etc., in both the text 
and the numbering of illustrations. In the text, ‘Eq. (1),” or 
“Eqs. (1) and (2)” are preferable to ‘‘Equation (1).”’ In cap- 
tions and legends, except for ‘‘Fig.’”’ and ‘‘Eq.”’ and table head- 
ings, write all words in full; do not abbreviate. Avoid placing 
explanatory written matter in the drawings; it should be in 
the text. 


MATHEMATICAL WorRK: Only the very simplest formulae 
should be typewritten; all others should be very carefully written 
in pen and ink, the writing to be large enough so that ample room 
is provided to mark mathematical matter for the printer. A 
considerable space for marking should be allowed above and be- 
low all equations. All complicated equations should be repeated 
on separate sheets with plenty of space left for marking. The 
solidus should be used for simple fractions appearing within 
the text. Make all expressions clear to the typesetter. Greek 
letters used in formulae should be clearly designated by name 
on the margin of the manuscript. All symbols should be clearly 
written and carefully checked. The difference between capital 
and lower-case letters should be clearly distinguished and care 
taken to avoid confusion between zero (0) and the letter (0), 
between the numeral (one) and the letter (ell) and the prime 
(‘), between alpha and a, kappa and k, u and mu, v and nu, n 
and eta. All subscripts and exponents should be clearly marked 
and.dots and bars over letters or mathematical expressions should 
be avoided. Avoid complicated exponents and subscripts. 
When it is necessary to repeat a complicated expression, it should 
be represented by some convenient symbol. 


NOMENCLATURE AND ABBREVIATIONS: The National Advisory 
Committee for Aeronautics Nomenclature should be used in pref- 
erence to any others. Standard abbreviations should be used, 
and it should be noted that most abbreviations are lower case, 
such as m.p.h., b.m.e.p., i.hp., b.hp., hp., ... ete. 
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